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Abstract 


Significant increases in sedimentation rate accompany the evolution of multicellularity. 
These increases should lead to rapid changes in ecological distribution, thereby affecting 
the costs and benefits of multicellularity and its likelihood to evolve. However, how genetic 
and cellular traits control this process, their likelihood of emergence over evolutionary time- 
scales, and the variation in these traits as multicellularity evolves are still poorly understood. 
Here, using isolates of the ichthyosporean genus Sphaeroforma-close unicellular relatives 
of animals with brief transient multicellular life stages-we demonstrate that sedimentation 
rate is a highly variable and evolvable trait affected by at least 2 distinct physical mecha- 
nisms. First, we find extensive (>300*) variation in sedimentation rates for different Sphaer- 
oforma species, mainly driven by size and density during the unicellular-to-multicellular life 
cycle transition. Second, using experimental evolution with sedimentation rate as a focal 
trait, we readily obtained, for the first time, fast settling and multicellular Sohaeroforma arc- 
tica isolates. Quantitative microscopy showed that increased sedimentation rates most 
often arose by incomplete cellular separation after cell division, leading to clonal “clumping” 
multicellular variants with increased size and density. Strikingly, density increases also 
arose by an acceleration of the nuclear doubling time relative to cell size. Similar size- and 
density-affecting phenotypes were observed in 4 additional species from the Sphaeroforma 
genus, suggesting that variation in these traits might be widespread in the marine habitat. 
By resequencing evolved isolates to high genomic coverage, we identified mutations in reg- 
ulators of cytokinesis, plasma membrane remodeling, and chromatin condensation that 
may contribute to both clump formation and the increase in the nuclear number-to-volume 
ratio. Taken together, this study illustrates how extensive cellular control of density and size 
drive sedimentation rate variation, likely shaping the onset and further evolution of 
multicellularity. 
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Introduction 


The emergence of multicellularity from single-celled life represents a major transition, which has 
occurred many times independently across the tree of life [1-8]. Multicellularity can arise either 
by aggregation of single cells that come together or from single cells that are maintained together 
clonally after division [9-11]. The unicellular and intermediate multicellular ancestors, which led 
to present-day multicellular organisms, have long been extinct [3], obscuring direct investigation 
of how multicellular life has emerged. However, several strategies have been used to study the 
emergence of multicellularity, including the use of experimental evolution (EE) approaches and 
the investigation of novel non-model organisms at pivotal positions in the tree of life. 

For example, EE under controlled conditions allows selection for diverse phenotypes 
[12,13], including multicellularity [14-18]. Using EE, Ratcliff and colleagues repeatedly 
observed the evolution of a simple form of multicellularity in Saccharomyces cerevisiae and 
Chlamydomonas reinhardtii in response to gravitational selection [14-19]. Similarly, multicel- 
lularity emerged in yeast as a mechanism to improve the use of public goods [17]. In all these 
cases, cells form clumps by incomplete separation of daughters from mother cells, instead of 
by postmitotic aggregation [20-22]. 

Alternatively, non-model organisms with key evolutionary positions can be used to better 
understand the emergence of multicellularity. In particular, the study of unicellular holozoans 
(Fig 1A), the closest unicellular relatives of animals, revealed that these organisms contain a 
rich repertoire of genes required for cell adhesion, cell signaling, and transcriptional regulation 
and that each unicellular holozoan lineage uses a distinct developmental mode that includes 
transient multicellular forms [2,23-27]. For instance, the choanoflagellate Salpingoeca rosetta 
can form clonal multicellular colonies through serial cell division in response to a sulphonoli- 
pid of bacterial origin [28-30], whereas the filasterean Capsaspora owczarzaki can form multi- 
cellular structures by aggregation [31]. Ichthyosporeans display a coenocytic life cycle unique 
among unicellular holozoan lineages and pass through a short and transient clonal multicellu- 
lar life stage prior to the release of newborn cells [32-36]. Despite their pivotal phylogenetic 
position, their rich “animal” genetic toolkit and the capacity to undergo transient multicellu- 
larity, to date, unicellular holozoans lineages and EE, have never been combined. 

Beside the formation of new biological structures and increases in organismal size, the 
emergence of multicellularity is frequently accompanied with an increase in sedimentation 
rate. Indeed, S. cerevisiae snowflake yeast, multicellular C. reinhardtii, and S. rosetta colonies 
sediment faster when compared to their unicellular counterparts [15,16]. Such correlation has 
been described in several marine phytoplankton species where multicellular life stages show 
faster sedimentation than unicellular ones [38-41]. This phenotype has a large impact on 
where in the water column these microbes proliferate [38,42-44] and thus is presumed to be 
under strong genetic control and selective pressure. Despite its capacity to affect the depth at 
which marine species flourish, the role of sedimentation rate, the potential impact of its varia- 
tion, and its connection to the emergence of multicellularity has not been systematically ana- 
lyzed across unicellular marine organisms, including in the closest unicellular relatives of 
animals with transient multicellular life stages. Here, using EE, we characterize how regulation 
of sedimentation rate can influence the emergence of stable multicellular life forms in the 
ichthyosporean Sphaeroforma genus, a close unicellular relative of animals. 


Results 


Sphaeroforma species exhibit large variation in sedimentation rates 


Similar to other ichthyosporeans, Sphaeroforma species proliferate through continuous rounds 
of nuclear divisions without cytokinesis to form a multinucleated coenocyte [34,35,45,46]. 
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Fig 1. Sedimentation dynamics of Sphaeroforma arctica coenocytes during the life cycle. (A) Cladogram representing the position of ichthyosporeans 
including Sphaeroforma species within the eukaryotic tree. (B) Schematic representation of the coenocytic life cycle of S. arctica. (C) Data derived from Fig 1 
[37] (gray points) were used to scale velocity measurements determined in our sedimentation assay to physical (m/s) units (red points) (Methods). Error bars 
represent the 95% CI for each unique genotype, time point, and temperature measurement presented in our study (N = 3 for each of 1 or 2 independent 
replications). Values were log-transformed prior to calculation of error. This figure is to illustrate where our data fit in the scheme of known plankton 
sedimentation rates. For our best estimations of cellular density and velocity in meters per second, a subset of this data from Smayda’s Appendix Table 1 was 
used (Methods). (D) Relying on the Smayda dataset as reference, we used measurements of sedimentation rate from our assay along with cellular perimeter 
measurements to calculate maximum likelihood estimates of excess cellular density. These estimates are plotted on a landscape illustrating the relationship 
between density and size on sedimentation rate (gray contour lines). Gray contour lines represent the predicted settling velocity of each genotype in pure water 
(excess density = 1,000 kg/m’) in units m/s. (E) Sedimentation rates of Sphaeroforma during the life cycle at 17°C. Every trace represents an independent 
experiment. (F) Distributions of nuclear content of S. arctica cells during the life cycle at 17°C, measured by microscopy (n > 500 cells per time point). (G) 
Quantification of mean DNA content per time point (measured by geometric mean) for cells grown in marine broth (MB) at 17°C (n > 500 cells per time 
point). (H) Average coenocyte volume per time point at 17°C (n = 100 coenocytes per time point). (I) Sedimentation rates of S. arctica coenocytes during the 
life cycle at 17°C. The data underlying this figure may be found at https://doi.org/10.6084/m9.figshare.18319232.v1. 


https://doi.org/10.1371/journal.pbio.3001551.g001 





Sphaeroforma coenocytes then undergo a coordinated cellularization process leading to the 
formation of a transient multicellular life stage resembling an epithelium [32]. This layer of 
cells then detaches and cell-walled newborn cells are released to the environment (Fig 1B) 
[32]. The entire life cycle prior to cellularization occurs in highly spherical (multinucleated) 
cells. 

Extensive literature has documented a positive correlation between cell size and sedimenta- 
tion rate, including during the life cycle of marine phytoplankton (Fig 1C) [37-41].This is con- 
sistent with Stoke’s law, which shows that the relationship between a spherical particle’s 
terminal sedimentation rate v in a fluid and its radius R should be determined by 


v=cp,R 


where c represents the scaled ratio of gravitational to viscosity constants, and py is the differ- 
ence between the particle’s and the fluid’s densities (see Methods). Therefore, even small shifts 
in a particle’s radius will lead to pronounced (i.e., quadratic) changes in sedimentation rate. 
Similarly, for particles sitting near the buoyancy threshold, small changes in density can lead 
to proportionally large changes in settling rate (Fig 1D). 
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Due to the nature of the coenocytic life cycle (Fig 1B), which is associated with an increase 
in the number of nuclei and coenocyte volume, we expected to observe an increase in cellular 
sedimentation rates over time [44,47,48].To better understand this relationship, throughout 
this study, we conducted overlapping experiments characterizing cell volume and sedimenta- 
tion rates of Sphaeroforma species for cultures growing at 17°C for 72 hours. For certain repli- 
cates of this core dataset, we included measurements of various genetic variants, the 
temperature dependence of phenotypes, fitness, as well as the speed of nuclear duplication. 
Overall, the measurements we report are highly reproducible with >95% variance explained 
for replicate measurements across phenotypes (S1A Fig and S1-S4 Data) (see Methods). This 
reflected a high heritability of the different phenotypes in a given environment. 

To begin, we measured the sedimentation rate of 5 different Sphaeroforma species that have 
been isolated from different habitats either free-living or host-derived, namely S. arctica, S. 
sirkka, S. napiecek, S. gastrica, and S. nootakensis, using a quantitative sedimentation rate assay 
based on changes in optical density over time [32,45,49-52]. According to these estimates, cell 
sedimentation rates varied greatly, from between 0.4 to 125 um per second (ie., up to 0.45 
meters per day) (Fig 1E). This broad variability over the life cycle and among Sphaeroforma 
species suggested that appreciable changes in size and/or cellular density should accompany 
different stages of the life cycle. 


The transient multicellular life stage of S. arctica is associated with an 
increase in sedimentation rate 


To better understand the cellular basis of sedimentation rate variation, we focused on S. arc- 
tica, the most studied Sphaeroforma species to date [32,34,50]. Using fixed cell imaging, we 
observed that synchronized cultures of S. arctica undergo a complete life cycle in about 48 
hours and can reach up to 256 nuclei per coenocyte before undergoing cellularization and 
releasing newborn cells (Fig 1F). Prior to this release, all cellular division occurs in highly 
spherical mother coenocytes. Consistent with previous results [34], nuclear division cycles 
were periodic during coenocytic growth and occurred on average every 9 to 10 hours, as mea- 
sured indirectly from changes in DNA content (Fig 1G). Average cell volume increased 
throughout the coenocytic cycle to reach its maximum value at 36 hours prior to the release of 
newborn cells (Fig 1H). Similarly, the sedimentation rate increased up to >300-fold the initial 
value after 36 hours (Fig 11). Altogether, we observe that every cycle of nuclear division is asso- 
ciated with a significant increase in nuclear content, volume, and sedimentation rate with a 
distinct peak just prior to cell release. As the transient multicellular life stage of S. arctica 
occurs during the latest stage of cellularization and ahead of cell release [32], our results sug- 
gest that it is tightly associated with an increased sedimentation rate. 


Experimental evolution of fast-settling mutants 


Given how cell size and division across the cell cycle are regulated, we reasoned that these vari- 
able traits should also be heritable and, hence, evolvable. To test this, we conducted an evolu- 
tion experiment to generate mutants with increased sedimentation rates [15,16,18]. Briefly, 10 
independent cultures of S. arctica (S1 to $10) originating from the same ancestral clone (AN) 
(Fig 2A) were diluted in fresh MB medium at 17°C (Fig 1C). Selection was performed every 24 
hours by allowing the cultures to sediment in tubes for 2 minutes before transferring and prop- 
agating the fastest-settling cells in fresh medium (Fig 2A). EE was continued for 56 transfers (8 
weeks) or about 364 nuclear generations, and a frozen stock was conserved every 7 transfers (1 
week) (see Methods) (Fig 2A, S1 Table). Most lineages exhibited a distinct clumping pheno- 
type to varying degrees across all evolved populations (Fig 2B). This phenotype has never been 
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Fig 2. Rapid evolution of fast-settling and clumpy mutants. (A) Schematic representation of the EE design using sedimentation as a selective pressure. All 
10 replicate populations (S1-S10) were isolated from the same ancestral culture (AN) and were subject to 56 transfers over 8 weeks of selection. (B) 
Representative images of the ancestral culture (AN) and 3 evolved mutants (S1, $5, and S6) with varying clumping capability after 5 days of growth. Bar, 

50 um. (C) Sedimentation rates of S. arctica evolved populations per number of transfers shows rapid emergence of fast-settling phenotypes, particularly in 
S1. (D) Sedimentation rates of single clones of S. arctica evolved mutants after 56 rounds of selection displays vast variations in sedimentation rate. 
Associated with S1 Video. (E) Head-to-head competitions between the clumpy evolved isolates S01 and S03 with the ancestral strain (AN) show that 
sedimentation rate is adaptive in the conditions of our evolution experiment. The rate of change of the frequency of the clumpy phenotype per day in log 2 
units was determined by least-squares regression. Error bars represent the 95% confidence interval based on 500 calculations of this statistic obtained by 
resamplings of the observed binomial count frequencies for each data point. The data underlying this graph may be found in S4 Data; see Methods. (F) 
Clump size distribution of evolved clones shows that fast-settling mutants (S1, $4, and S9) show bigger clump size (n > 800 measurements per time point). 
The data underlying this figure may be found at https://doi.org/10.6084/m9.figshare.18319232.v1. EE, experimental evolution. 





https://doi.org/10.1371/journal.pbio.3001551.g002 


observed in any Sphaeroforma culture, despite variation of culturing conditions across labs, 
and the hundreds of routine passages without selection, thus suggesting that the appearance of 
such phenotype is linked to our selection regime. The clumps of cells were not maintained by 
ionic forces or protein-dependent interactions, however could be separated by mild sonication 
without leading to cell lysis (S2B Fig). To assess when sedimentation rates increased during the 
selection process for each evolved population, we synchronized cultures using mild sonication 
before dilution in fresh medium and allowed them to undergo a complete life cycle before 
measuring the sedimentation rate (Fig 2C). We observed that populations $1, $4, and S9 had 
the highest sedimentation rates at the end of the evolution experiment (Fig 2C). We observed 
a dramatic increase in sedimentation rate already after 14 transfers for population S1 (approxi- 
mately 91 generations) (Fig 2C). To assess and compare the variability in sedimentation rates 
among evolved lineages, we isolated and characterized a single clone from each evolved cul- 
ture. Our results show that evolved clones all settled significantly faster than the common 
ancestor but that there was stark variation in sedimentation rates at the end of the EE (Fig 2D). 
Isolates from lineages $1, S4, and S9 settled the fastest upon sedimentation (S1 Video), while 
clones from lineages S2, S3, S5, and S7 were intermediate, and lineages S6, S8, and S10 settled 
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the slowest (Fig 2D). Taken together, our results show that we can rapidly and routinely evolve 
fast-settling mutants in S. arctica using EE but that the outcome is not uniform across lineages. 

To assess whether fast sedimentation is adaptive in our EE setup, we independently com- 
peted 2 fast-sedimenting isolates, S01 and S03, against their common ancestor (AN). For this 
relative fitness assay, we mixed monocultures of either S01 or S03 with AN at an approximately 
1:1 ratio and subjected them to the same selective regime as during EE. Parallel cultures of 
each 1:1 mix grown and maintained without sedimentation selection served as experimental 
control (see Methods). Our results show that the proportion of clumpy S01 and S03 cells both 
significantly increased over time (a log2 selection coefficient of 0.23 and 0.13 per day for S01 
and S03, respectively, p < le-7), suggesting that this trait was adaptive in our experimental 
setup (Figs 2E, S1A, and S2D; S4 Data). In contrast, S01 and S03 cells decreased in frequency 
relative to the AN strain in the control growth experiment without sedimentation (a log2 selec- 
tion coefficient of —0.13 and —0.15 per day for S01 and S03, respectively, p < le-7), suggesting 
that the heritable changes that these strains bore were deleterious in normal laboratory propa- 
gation conditions (Figs 2E and S2D). 

We next sought to characterize the cellular mechanisms giving rise to variation in cellular 
sedimentation rates in evolved clones. Using quantitative microscopy, we show that clump 
perimeter and average number of cells observed per clump correlated highly with sedimenta- 
tion rates (Figs 2E, S2C, and S2E). Indeed, across all experiments reported in this study, we 
found that using Stokes’ law, a single globally fitted density parameter and our size measure- 
ments, could explain the majority of variance in sedimentation rate (R” = 0.69, RMSE (as a 
proportion of the range in observations) = 18%). 


Fast-sedimenting mutants form clonal clumps 


Earlier, we defined 3 distinct developmental stages of S. arctica life cycle: (i) cell growth with 
an increase in coenocyte volume; (ii) cellularization, which coincides with actomyosin network 
formation and plasma membrane invaginations; and (iii) release of newborn cells (Fig 1B) 
[32]. A key developmental movement named “flip,” defined by an abrupt internal morphologi- 
cal change in the coenocyte, can be used as a reference point to characterize life cycle stages. 
Prior to this event (pre-flip), actomyosin-dependent plasma membrane invaginations occur, 
while afterwards (post-flip), the cell wall is formed around individual cells prior to their release 
from the coenocyte [32]. Using time-lapse microscopy at 12°C, we observed that clump forma- 
tion in fast-settling mutants coincides with cell release (Fig 3A, S2 Video). 

Importantly, no new cell aggregation processes were detected after cell release, suggesting 
that all clumps formed either prior to or concomitantly with the release of newborn cells (S2 
Video). To examine whether clumps formed due to defects at the level of either plasma mem- 
brane or cell wall, we stained plasma membranes using FM4-64 and cell walls using calcofluor- 
white. We found that the process of plasma membrane invaginations during cellularization 
appears to be unchanged prior to flip (Fig 3B, S3 Video) and that all newborn cells, even in the 
clumps, have a distinctive cell wall surrounding them (Figs 3C and S3A). These results show 
that clumps are formed post-flip in fast-settling mutants. 

As no cells appear to aggregate after release of newborn cells, our results suggest that clumps 
are maintained together in a clonal form. To further confirm this, we sonicated both the ances- 
tor (AN) and S1 clumps and stained them separately with 2 distinct fluorescent dyes prior to 
mixing them in a 1:1 ratio (S3B Fig). After 1 complete life cycle, we observed that out of 198 
clumps, 183 contained only evolved S1 derivates, whereas ancestral cells only randomly and 
sporadically formed clumps under the conditions of the experiment (a total of 12 small 
clumps). Only 3 of the clumps (approximately 1.5%) contained cells of both colors (S3B Fig), 
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Fig 3. Clumps are formed by incomplete cell-cell separation in fast-settling mutants. (A) Time-lapse images of the 
life cycle of S. arctica ancestral strain (AN) and the fast-settling mutants (S1, S4, and S9) at 12°C show that clumps are 
observed concomitantly with release of newborn cells. Associated with $2 Video. Bar, 50 um. (B) Plasma membrane 
staining using FM4-64 show that plasma membrane invaginations during cellularization seems to occur normally in 
fast-settling mutants. Associated with $3 Video. Bar, 50 um. (C) Cell wall staining using calcofluor-white indicates that 
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cells are still separated by a cell wall inside fast-settling clumps, suggesting that clumps are formed post-flip. Bar, 

50 um. (D) Number of cells detaching per clumps at cell release and for the following 3 hours, measured from time- 
lapse movies (n = 50 coenocytes at cell release per strain). (E) Coenocyte volume at flip, measured from time-lapse 
movies, shows that S4 and S9 are significantly smaller when compared to AN and S1 (n = 50 coenocytes at flip per 
strain). The data underlying this figure may be found at https://doi.org/10.6084/m9.figshare.18319232.v1. 


https://doi.org/10.1371/journal.pbio.3001551.g003 


and each only contained a single AN cell trapped inside a smaller $1 clump. In a control exper- 
iment with two differentially stained S1 cultures, we observed an almost identically small num- 
ber (i.e., approximately 1.7%) of mixed clumps (S3B Fig); hence, it appears that mixes of cells 
happen only sporadically at low frequencies by random association, irrespective of genotype. 
Altogether, our results suggest that the evolved clump phenotype is not a result of spontaneous 
cell aggregation and instead arises from incomplete detachment between cells. 

We next examined the kinetics of cellularization and the process by which cells propagate 
as part of clumps. We counted the number of cells detaching from clumps at cell release and 
observed variable detachment across all 3 fast-settling mutants. The clumpiest isolate, popula- 
tion S1, was the least prone to cell detachment (Fig 3A and 3D, S2 and S3 Videos). Intriguingly, 
despite exhibiting similar sedimentation rates, 4 and S9 clumps showed a higher detachment 
frequency than S1 (Fig 3A and 3D). Image analysis of time-lapse movies also showed that life 
stage durations varied among fast-settling mutants at 12°C, with S1, $4, and S9 initiating cellu- 
larization and undergoing flip, respectively, 5.5, 7.5, and 10.5 hours earlier than the ancestor 
(S3C and S3D Fig). Post-flip duration also varied significantly among mutants, with $1 and S9 
requiring more time to release cells compared to $4 (S3C and S3D Fig). While many aspects of 
the replication cycle dynamics were variable for these mutants, the duration of cellularization 
was fairly invariant (approximately 9 hours) (S3C and S3D Fig). Finally, measurements of coe- 
nocyte volume show that $4 and S9 coenocytes undergo flip at substantially smaller volumes 
(approximately 1.8 and approximately 3.3 times smaller, respectively) compared to AN or S1 
(Fig 3E). These results show that, despite their shared capacity to clump and similar sedimen- 
tation rates, fast-settling mutants exhibit significant variability in their life cycle dynamics, 
with S4 and S9 mutants initiating cellularization earlier, dispersing from clumps with a higher 
frequency, and undergoing flip and cell release at smaller coenocyte volumes compared to the 
S1 mutant. 


Increased nuclear number-to-volume ratio leads to faster sedimentation 


Above, we observed that $4 and S9 mutants can sediment as fast as S1 despite their smaller 
coenocyte volumes, suggesting an alternative regulation mechanism of sedimentation rate. 
Across all experiments reported in this study, we found that approximately 31% of the variance 
in observations of sedimentation rate could not be explained by cell size only, suggesting that 
cellular density might also contribute to this variation [37,40,41]. From Stoke’s law, we calcu- 
lated that excess cellular density, i.e., cellular density minus that of distilled water (1,000 kg/ 
m°), might vary between 40 and 300 kg/m’ for S. arctica wild-type and evolved clones across 
their life cycle—the upper limits between the densities of pure protein and pure cellulose. Val- 
ues reached approximately 650 kg/m’ for wild S. nootakensis (Snoo) soon after cellularization, 
approaching the excess density of pure nucleic acid (Fig 1D). During cell cycle stages prior to 
cellularization, when cells were most spherical (<36 hours), excess density varied from 40 to 
200 kg/m? across Sphaeroforma isolates. 

To better characterize the relationship between sedimentation rate, cell cycle, and size, we 
performed higher resolution measurements of sedimentation rates over the complete life cycle 
of the ancestor and all 3 fast-settling mutants at 12°C and 17°C. Consistent with their capacity 
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to form clumps, we observed that the sedimentation rate of all fast-settling mutants increases 
during growth but, unlike the ancestor, does not recover after cell release to their original levels 
(Figs 4A and S4A). Interestingly, we noticed that individual $4 and S9 coenocytes sediment 
faster (approximately 2.5x and approximately 1.6x, respectively) than $1 or AN even before 
clump formation (24 to 36 hours time points) (Figs 4A and S4A). Such increase in sedimenta- 
tion rate was not due to a rise in cell size or change in cell shape as both $4 and S9 exhibit 
smaller cell perimeters throughout the cell cycle (Figs 4B, 4C, S4B, and S4C). Rather, excess 
cellular density estimations show that both $4 and S9, even prior to cell release and clump for- 
mation, tend to be on average 3x denser when compared to the ancestor (Figs 4D and S4D). 
Altogether, these results show that both cell size and cell density contribute to sedimentation 
rate variation in S. arctica. 

As cell size and nuclear division cycles are decoupled in S. arctica [34], we reasoned that 
increased cell density in $4 and S9 could be caused by an acceleration of nuclear divisions lead- 
ing to a rise in the number of nuclei per volume. Using DAPI staining to label nuclear DNA, 
we observed that $4 and S9 undergo nuclear duplication faster (approximately 2 hours) than 
both AN and S1 (Figs 4E, 4F, and S4E-S4H). By carefully examining the volumes of coeno- 
cytes containing the same number of nuclei at the single-cell level, we show that for the same 
nuclear content, $4 and S9 tend to be 30% to 45% smaller in volume when compared to the 
ancestor (Figs 4G and S4I). Consequently, both $4 and S9 exhibited the highest number of 
nuclei per volume (nuclear number-to-volume ratio) (Figs 4H and S4J). Taken together, these 
results argue that cell density can contribute an appreciable amount to cellular sedimentation 
rates (up to approximately 50 um/s) and that mechanistically this could arise by faster nuclear 
doubling times relative to cell size. 


Evolved genetic variation correlating with fast sedimentation 


Up to now, our results suggest that S. arctica mutants evolved faster sedimentation using 2 
strategies: (i) clump formation; and (ii) increased nuclear number-to-volume ratio. We found 
that sedimentation rate variation was highly heritable, persisting for 780 generations of passag- 
ing for all 10 isolates without selection for sedimentation phenotype, suggesting that the phe- 
notypes have a genetic basis. To test this, we resequenced the whole genomes of both the 
ancestral clone (AN) and 1 evolved clone per lineage (S1 to $10) obtained at the conclusion of 
the evolution experiment (Week 8) with high (>30-fold) coverage of the very large genome 
size of S. arctica, at 143 Mbp. Following very careful variant filtering (S2 and S3 Tables), we 
identified a total of only 26 independently evolved variants with an average of 2.6 mutations 
per clone (range 1 to 5 per clone) (S4 Table). Of the 26 variants, 24 (approximately 92.3%) 
were SNPs (11 coding, and 13 intergenic or intronic), and 2 were insertions (1 coding, 1 inter- 
genic) (S4 Table). Beside the fact that 2 of these variants are identical SNPs at the same position 
(Sarc4_g3900) and have independently evolved in 2 different backgrounds (S4 Table), no 
other form of genetic parallelism was detected. While it is very likely that some of these muta- 
tions hitchhiked in the background of beneficial driver mutations, we found that the coding 
SNP variants were clearly skewed toward nonsynonymous changes (9:2), with a cumulative 
dN:dS ratio of 1.32. This indicates the general presence of positive selection and, hence, adap- 
tative evolution driving this molecular pattern. Based on COG assignments (S6 Table), 4 of the 
changes are orthologous to genes implicated in signal transduction, and 2 are related to genes 
with functions in DNA binding or chromosome condensation. 

In the absence of molecular genetic tools and functional knowledge of many of the mutated 
gene targets, we set out to better understand how the distinct genetic variants could have influ- 
enced sedimentation rates and clump formation by examining the predicted expression 





PLOS Biology | https://doi.org/10.1371/journal.pbio.3001551 March 29, 2022 9/31 


PLOS BIOLOGY 


Regulation of sedimentation rate shapes the evolution of multicellularity in a close relative of animals 





























=1,4x10° 
A Cc D p=6,1x105 
300 250 200 
¢ a asixio* —pataxto” ANE J A 
8 —=S1 E € 580 PEO Ra hexit as1 p=1,7x104 
= Se = s p=4,3x10® —-P=4,5x107 BS4 5150 p=2,1x10* 9-4 3,198 
= 20089 2 © =0,5 a aso %@ 2 oe 
2 2 ® 150 as a 5 a a 
= = = Hl & 100 o | 
5 a a | & r f 
® 400 ® o 5 | 2 p=0,66 
2 2 fo) is] [ay —a BAN 
o 3 gS 2 ° 50-88 = 
£ > 2 50-5 6 5 2 7 s1 
3 <x x i] i: g 5 is o BS4 
| 
o _ 6 p=5,9x107 9 1 asg 
0 12 24 36 48 60 72 0 12 24 36 48 60 72 60 24 60 
Time (h) Time (h) Time (h) Time (h) 
E F G p=1,7x10” . Hes pssst bales 
32 — 10-9 5 _ 100 ——— , # a 4 6 H 
~ | € 2 21 p=8,9x10" 5 a 
< = p=2,2x10 . 
= © = = Xx p=7x10- ra 2 
eS 5 p08 - 10 —-—+ i ° us 4 p=7x104 ot p=0,1 ¢ 
=e £8 eo 2 yy a p=0,001 = FOr es 
62 7 a oo § © p-o001 TT ; : \ ; f 4 
85 3 2 — p=6,9x10"" p=2,extov Sk 8 8 
ge4 37 es > 2 i 
zé SS 2 4 il 8 ———, BAN 24; 
a8 7 p=9,6x10* S =1,2x10°7 
2 2 3 6 8 T A 7 i as1 
3 p=6,1x10° 8 - —S4 
if 5 S04 aso ly 
0 12 24 36 48 60 72 AN S1 S4 $9 1 4 16 64 ANS1S4 S9 
Time (h) DNA content 
120 | 31 200-| 54 140-| gg “30 
E 100. 5 5 
= | = 150 = 105 
8 804 : 2 
oO oO oO 
© 60 FE z 
5 60 5 100 5 70 15 
a Qa Re} 
oO oO oO 
gs 40- sg = 
Zz | Zz Zz 
x | ma 50 we 3 
E 20- 5 7 ~~ = 
Se, 
0 le 2 2 eo . 2 e 2° © o 0 eee _o_o__o_e_eoe 0 0 
a et ie a 7 
12 18 24 30 36 42 48 54 60 66 12 18 24 30 36 42 48 54 60 66 12 18 24 30 36 42 48 54 60 66 


Time (h) 


—e~ Sarc4_g12019!} Ist1 
—e~ Sarc4_g3900 


Time (h) 
—e— Sarc4_g11520 
Dehydrocholesterol reductase 
—e— Sarc4_g32431 |! SLK 


Time (h) 


—e— Sarc4_g3019! RCC1 
_,_ Sare4 | 97653* 
bHLH Transcription factor 


Growth —e- Sarc4_g23215 2 
m Pre-flip a3 — Saro4_g33270 eaves oT i880 a 
Post-fip Cellularization —— Saro4_g33124 C2H2 Zinc finger Transcription factor 


—e Sarc4_g22580* 


Cell rel 
H Cell release —e Sarc4_g11957* 


—e Sarc4_g11958* 


Fig 4. Sedimentation rates variation in fast-settling mutants is associated with variation in cell size and cellular density. (A) Sedimentation rates of S. 
arctica AN and evolved mutants during the life cycle at 17°C. Every trace represents an independent experiment. (B) Average perimeter measured from 
fixed cells every 12 hours over a complete life cycle of 72 hours at 17°C shows that fast-settling mutants increase their size upon cell release. Every trace 
represents an independent experiment (n > 180 measurements per time point for each independent experiment). (C) Average perimeter of fast-settling cells 
and clumps at 24 and 60 hours, respectively, show that $4 and S9 cells and clumps have a smaller size when compared to S1. Every square represents an 
independent experiment, and the white circle represents the median (n > 180 coenocytes per time point for each independent experiment). (D) Excess 
cellular density of fast-settling individual coenocytes (before cellularization) and clumps (after cell release) at 24 hours and 60 hours, respectively, show that 
S4 and S9 single cells are denser when compared to S1 and AN. Every square represents an independent experiment, and the white circle represents the 
median. (E) Quantification of mean DNA content per time point for fast-settling mutants grown in MB at 17°C. Every trace represents an independent 
experiment (1 > 400 coenocytes per time point for each independent experiment). (F) Nuclear doubling time, calculated by linear regression of mean 
nuclear content at time points from 0 hour to 24 hours. Every square represents an independent experiment, and the white circle represents the median 

(n > 400 coenocytes per time point for each independent experiment). (G) Boxplots of cell volume measurements of DAPI-stained fixed cells. For 1-, 4-, 16-, 
and 64-nuclei cells. Cells with 1 nucleus represent newborn cells at the end of the experiment (” > 80 coenocytes per DNA content). (H) Boxplots of nuclear 
number-to-volume ratio of DAPI-stained cells show significant increase for $4 and S9 fast-settling mutants. Every square represents an independent 
experiment, and the white circle represents the median (n > 600 coenocytes per strain). (I-K) Temporal transcript abundance of genes mutated in fast- 
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settling phenotype across the native life cycle of S. arctica. The data underlying these graphs may be found in S5 Table. The data underlying this figure may 
be found at https://doi.org/10.6084/m9.figshare.18319232.v1. 


https://doi.org/10.1371/journal.pbio.3001551.g004 


dynamics of mutated genes across the cell cycle. The data were derived from a recently pub- 
lished time-resolved transcriptomics dataset of the S. arctica life cycle [32]. Of the mutated 
genes, 9 showed no expression during the native life cycle, whereas 12 displayed dynamical 
expression during cellularization, and the remaining 5 genes were more or less stably 
expressed (Figs 4I-4K and S4K, S5 Table). We also annotated all mutation-associated genes 
based on a recent comprehensive orthology search [53] (S6 Table). 

The fastest-settling and clumpiest mutant isolated from the population S1 bore a synonymous 
mutation of a homolog of increased sodium tolerance 1 (Ist1) superfamily and as such likely 
impacts gene expression rather than gene function. This gene shows a dynamic expression during 
cellularization and codes for a conserved protein involved in multivesicular body (MVB) protein 
sorting (Fig 41, S6 Table) [54,55]. In humans, hIST1, also known as KIAA0174, is a regulator of 
the endosomal sorting complex required for transport (ESCRT) pathway and has been shown to 
be essential for cytokinesis in mammalian cells [56]. Similarly, Ist1 orthologs in both budding and 
fission yeasts play a role in MVB sorting pathway and, when deleted, exhibit a multiseptated phe- 
notype consistent with a role in cytokinesis and cell separation [57,58]. 

In the clone derived from population S4, we observed 5 distinct mutations, 2 nonsynon- 
ymous, 1 intergenic, and 2 intronic SNPs. Among the 2 nonsynonymous SNPs, one causes a 
E90G change in a homolog of human Kanadaptin (SLC4A1AP), which may play a role in sig- 
nal transduction [59,60]. Among the noncoding SNPs, 1 mutation is found in an intron of the 
7-dehydrocholesterol reductase (DHCR7), expressed during cellularization and known to be 
key in the cholesterol biosynthesis pathway [61,62], and the second intronic mutation codes 
for a STE20-like kinase (SLK), which plays numerous roles in cell cycle signaling and actin 
cytoskeleton regulation (Fig 4J) [63-66]. 

Finally, among the 5 mutations discovered in the clone derived from population S9, 2 muta- 
tions are in transcription factors that are continually expressed during the cell cycle: an intronic 
SNP in a basic helix-loop-helix (bHLH) transcription factor, and the sole nonsynonymous 
SNP leading to A923V change in a gene predicted to encode a nucleotide binding C2H2 Zn fin- 
ger domain [67]. A third mutation was found in an intron of the highly and dynamically 
expressed homolog of the regulator of chromosome condensation 1 (RCC1; Fig 4K) [68-70]. 
RCC1 is a chromatin-associated protein implicated in several processes including nuclear for- 
mation, mRNA splicing, and DNA replication [69,71-74]. Thus, it may contribute to the accel- 
erated nuclear duplication cycle observed in S9 by impacting cell cycle progression. Altogether, 
the mutations identified in both $4 and S9 may affect both cellularization and cell separation. 

Among the variants detected in evolved clones with intermediate-settling phenotype, we 
highlight an intergenic mutation 118 bp downstream of Dynamin-1 known to be essential for 
cytokinesis across different taxa [75-77], and a nonsynonymous mutation in a protein similar 
to Fibrillin-2 (Sarc4_g7365T), which is an extracellular matrix (ECM) glycoprotein essential for 
the formation of elastic fibers in animals (S4K Fig) [78-80]. Altogether, our results across all iso- 
lates suggest that a large mutational target affects cellular sedimentation and multicellularity. 


Sedimentation rate variation across Sphaeroforma species is driven by cell 
size and density 


Lastly, we examined whether the variation in cell sedimentation observed across different 
Sphaeroforma species (Fig 1E) could also be explained by clumping or increased nuclear 
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number-to-volume ratio. To do so, we investigated the life cycle dynamics, coenocyte volume, 
and nuclear duplication time among distinct Sphaeroforma sister species. To date, 6 different 
Sphaeroforma species have been isolated either in a free-living form or derived from different 
marine hosts: S. arctica, S. sirkka, S. napiecek, S. tapetis, S. gastrica, and S. nootakensis (Fig 5A) 
[32,46,49-52]. Using previously established growth methods for S. arctica combined with live 
and fixed imaging, we first observed that all sister species but S. tapetis show a synchronized 
coenocytic life cycle (Figs 5B and S5A, S4 Video). Similar to abovementioned results with fast- 
settling mutants, sedimentation rate variations (Fig 1E) could be explained by variations in 
both cell size and cellular density. Indeed, we first observed that both S. gastrica and S. noota- 
kensis occasionally form clumps, exhibit a lower frequency of cell detachment, and thus have 
an increased cellular density after cell release compared to the other Sphaeroforma species 
(Figs 1C, 5B, and 5C, S4 Video). Additionally, despite their increased sedimentation rate (Fig 
1C), we found that all sister species exhibited approximately 20% to 45% smaller coenocyte 
size prior to cell release when compared to S. arctica, which reveals an increased cellular den- 
sity (Figs 5D, S5B, and S5D, S4 Video). Similar to the fast-settling mutants S4 and S9 above, 
the increase in cellular density was associated with an acceleration of the nuclear division 
cycles and the subsequent rise in the nuclear number-to-volume ratio (Figs 5D-5F and S5D- 
S5F). Notably, S. sirkka and S. napiecek, both previously isolated as free-living [51,52], exhibit 
an increase in nuclear number-to-volume ratio but no ability to form clumps. Altogether, our 
results show that, similarly to experimentally evolved strains, fast sedimentation variation 
could occur by both clump formation and/or increase in the nuclear number-to-volume ratio 
for Sphaeroforma species (Fig 5H) and thus might represent in itself a widespread and highly 
variable phenotype in the marine habitat. 


Discussion 


Here, we performed, for the first time, EE on one of the closest unicellular relatives of animals, 
and we demonstrate that, under suitable selection pressure, the ichthyosporean S. arctica can 
evolve stable multicellularity. In particular, we observed the independent rise of clump forma- 
tion, and faster settling phenotypes across populations within less than 400 generations. The 
precise detectable onset of the phenotypes varied across lineages and occurred as early as 
approximately 91 generations in lineage S1. We further show that faster sedimentation pheno- 
types are highly heritable and, according to direct competition assays with the common ances- 
tor, are expected to be adaptive under the environmental conditions of the experiment. Our 
results add to previous observations of the rapid emergence of multicellularity in yeast and 
green algae, which all can evolve multicellular clump-forming structures within short evolu- 
tionary timescales [15,16,18]. As ichthyosporeans proliferate through an uncommon coeno- 
cytic life cycle [32,33,35,45,46,51], our results show that this evolutionary process of selection 
for faster sedimentation is accessible at microevolutionary timescales across taxa and organ- 
isms with highly diverged modes of proliferation. 

In this study, we show that all fast-settling S. arctica cells increased their cell size by increas- 
ing cell adhesion post-cellularization, leading to the formation of clumps (Fig 5H). Such results 
are analogous to cell cluster formation in snowflake yeast and C. reinhardtii, which arises 
through incomplete separation of mother and daughter cells [15-17]. Altogether, these results 
suggest that regulation of sedimentation rate can constrain unicellular species to generate mul- 
ticellular cell phenotypes by increasing their cell adhesion efficiency. However, we found that 
approximately 31% of the variance could not be explained by cell (clump) size only. Indeed, 2 
fast-settling mutants (S4 and S9) exhibited an increase in sedimentation rate prior to clump 
formation, which was associated with an accelerated nuclear division cell cycle leading to an 
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Fig 5. Sedimentation rate variation across Sphaeroforma species is associated with clumping and increased nuclear number-to-volume ratio. (A) A 
cladogram [49], representing the position of all Sphaeroforma sister species used in the study and their observed phenotypes relative to S. arctica. Cells exhibit 
either variation in clumpiness and/or nuclear number-to-volume ratio (N/C), which is regulated by the doubling time relative to volume. (B) Representative 
images of different Sphaeroforma sp. at cell release. Arrowheads indicate the formation of clumps. Associated with S4 Video. Bar, 50 um. (C) Number of cells 
detaching at cell release, measured from time-lapse movies, show significant differences among the different sister species. (n > 48 coenocytes per 
Sphaeroforma sp.). (D) Excess cellular density of Sphaeroforma sp. cells (before cellularization) and clumps (after cell release) at 36 and 72 hours, respectively. 
Every square represents an independent experiment, and the white circle represents the median. (E) Quantification of mean DNA content per time point for 
Sphaeroforma sp. grown in MB at 17°C. Every trace represents an independent experiment (n > 300 coenocytes/time point for each independent experiment). 
(F) Nuclear doubling time, calculated by linear regression of mean nuclear content at time points from 0 hour to 24 hours at 17°C. Every square represents an 
independent experiment, and the white circle represents the median. (n > 300 coenocytes per time point for each independent experiment). (G) Boxplots of 
nuclear number-to-volume ratio of DAPI-stained cells at 17°C for Sphaeroforma sister species. Every square represents an independent experiment, and the 
white circle represents the median. (n > 300 coenocytes per strain). (H) Scheme illustrating a landscape of how phenotypic and genotypic plasticity varies over 
the physical dimensions of size and cellular density. Gray axis shows the cellular traits affecting these physical dimensions. Colors represent how the different 
Sphaeroforma species investigated in this study behave across their life cycle. Black = S. arctica (AN), Blue = S. arctica (S1 mutant), Pink = S. arctica (S4 and S9 
mutants, S. gastrica, S. nootakensis), Green = S. sirkka, S. napiecek. Thick lines represent the typical cellular behavior. Thin and dotted lines indicate that 
clumping cells may restart the process at any position. The data underlying this figure may be found at https://doi.org/10.6084/m9.figshare.18319232.v1. 


https://doi.org/10.1371/journal.pbio.3001551.g005 
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increase in the number of nuclei per unit volume (Figs 4G and 5H). Previous results have 
shown that, in S. arctica, both nuclear duplication cycles and cell size are uncoupled [34]. Our 
results support these findings and indicate that nuclear division cycles and cell size could be 
regulated separately, allowing adaptive change in either, and independently of one another. 

By analyzing the genomes of evolved isolates, we identified an overabundance of nonsynon- 
ymous mutations, indicating positive selection, with at least 1 and up to 5 independent muta- 
tions in each lineage. Many of the better-characterized genes that carry mutations in either 
coding or intergenic regions are dynamically expressed during cellularization. Several muta- 
tions were found in genes coding for cytoskeletal regulators and cytokinesis proteins, which is 
consistent with previous studies in which cytokinesis-deficient mutants were often associated 
with an incomplete cell separation across taxa [81-85]. Other mutations were found in genes 
involved in cell signaling, plasma membrane remodeling, and chromatin condensation regula- 
tors reflecting a large and accessible mutational target affecting sedimentation rate phenotypes. 
Cytokinesis defects may lead to incomplete cell separation, which may explain part of the cell 
clumping phenotype observed in the fast sedimentors. Given that similar phenotypes emerged 
independently multiple times with no genetic overlap, we conclude that the mutational target 
for these traits could be quite large for S. arctica, opening the possibility for variation and evo- 
lution in multicellularity-related phenotypes. 

We found that closely related species exhibit widespread variation in both clump formation 
and nuclei-to-volume ratio. Given the evolvability of the trait and the few natural samples, we 
examined it is difficult to argue which cellular states affecting sedimentation rates are ancestral. 
For example, considering the species tree and the trait of clumpiness (Fig 5A and 5B, S4 
Video), only 3 Sphaerforma species exhibited this behavior under our experimental conditions: 
S. gastrica, S. tapetis, and S. nootkatensis. Despite characterizing for the first time all the avail- 
able isolates of Sphaeroforma sp., the number of studied isolates remains limited (6 species), 
and, thus, we cannot distinguish whether this trait is “ancestral” or “novel.” Rather, our experi- 
mental data suggest that clumpiness and density could evolve over short microevolutionary 
timescales such as those we have measured in the lab and this in both directions. Evolved lines 
with highly clumping phenotypes had the competitive edge over their nonclumping common 
ancestor in relative fitness assays. In contrast, during exponential growth and in the absence of 
selection, evolved clumpy line were selectively disfavoured relative to their common ancestor, 
which suggests that clumpy phenotypes are prone to rapid replacement by less clumpy or 
dense morphs when environments change. 

Indeed, marine organisms have evolved various passive or active means of maintaining 
their position in the water column, for example, using motility and/or ingenious approaches 
to regulate buoyancy [44,86-90]. Sphaeroforma species are remarkably spherical, immobile, 
lack flagella, and yet exhibit a substantial increase in cell size and density over the life cycle, 
thus representing a challenge to maintaining buoyancy in marine habitat. This work estab- 
lishes that Sphaeroforma’s cell size and density are subject to tight cellular control and are 
highly evolvable traits. Taken together, these observations suggest that sedimentation rate is a 
highly evolvable trait, which itself likely shapes the gain and loss of multicellularity. 


Methods 
Culture conditions 


All Sphaeroforma sp. cultures were grown and synchronized as described previously for S. arc- 
tica [32,34]. Briefly, saturated cultures in MB (Difco BD, NJ, USA; 37.4 g/L) were diluted into 
fresh medium at low density (1:250 dilution of the saturated culture) and grown in rectangular 
canted neck cell culture flask with vented cap (Falcon; ref: 353108) at either 17°C or 12°C, 
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resulting in a synchronously growing culture. Saturated culture of Sphaeroforma sp. are 
obtained after 3 weeks of growth in MB. 


Experimental evolution 


Ten replicate population (S1 to S10) of genetically identical S. arctica (AN) were first diluted 
250-fold in 5 ml of MB and grown in rectangular canted neck cell culture flask with vented cap 
(Falcon; ref: 353108) at 17°C. Cells were grown at 17°C rather than the previously used 12°C in 
order to increase growth rates and accelerate evolutionary outcome. Every 24 hours, the entire 
population was transferred into a 15-ml falcon tube and allowed to sediment for 2 minutes on the 
bench. A 5-ml pipette was then positioned vertically and used to collect 500 ul of cell culture from 
the bottom of the falcon tube. The cells, still vertically positioned in the pipette, were then allowed 
to sediment once more for 15 seconds before the transfer of a single drop, equivalent to 20 ul, into 
5 ml of fresh MB (approximately 250x dilution). Every 7 transfers, a frozen fossil was conserved 
by adding 10% of DMSO to 1 ml of culture and preserved at —80°C. Single clones of each replicate 
population (S1 to $10) were obtained at the end of week 8 by serial dilutions. Since S. arctica 
grows as a coenocyte, a temporal generation is not defined by a complete coenocytic cycle but is 
equivalent to the doubling of the number of nuclei. We estimated the nuclear doubling time by 
measuring the number of cells and the number of nuclei per cell for each transfer separately (S1 
Table). Briefly, the entire EE experiment comprised approximately 364 generations, in which all 
populations underwent a total of 28 complete coenocytic cycles. Each coenocytic cycle included 
two 24 hours subpassages and comprised a total of approximately 13 doublings (S2A Fig, S1 
Table). Importantly, the effective population size was kept in a very narrow range across subpas- 
sages (at approximately 10°) and thus over the entire experiment (S2A Fig, S1 Table). Therefore, 
in evolutionary terms, the population size was consistently high enough to favor of natural selec- 
tion over random evolution throughout the course of the experiment. This assumption was recon- 
firmed genetically by deriving a dN:dS ratio >1 from sequencing data (see main text). 


Sedimentation rate measurements 


Sedimentation rate was measured for Sphaeroforma sp. every 12 hours for a total of 72 hours 
unless indicated otherwise. To ensure reproducibility and homogeneous results, saturated 
Sphaeroforma sp. cultures were sonicated prior to the dilution in fresh MB media (250-fold 
dilution) using a Branson 450 Digital Sonifier (3 pulses of 15 seconds, 10% amplitude). For 
each measurement, obtained either from different stages of the cell cycle or from different 
Sphaeroforma species, 1 ml of cell culture was added into a disposable plastic spectrophotome- 
ter cuvette (semi-micro, 1.5 ml) and homogenized by vortex. Optical density (OD¢o0) was 
measured using an Eppendorf Biophotometer (Model #1631) at T = 0, corresponding to the 
first time point after placing the cuvette in the spectrophotometer. The OD¢o9 was then contin- 
uously measured every 30 seconds for 3 minutes while cells were slowly sedimenting in the 
cuvette. To ensure that OD¢o9 measurements stayed within the detection limits of the spectro- 
photometer, early life stages (TO-T48) were not diluted in the cuvette, whereas later life stages 
(T60-T72) were diluted 1/100 in fresh MB media. 

For assessing clump dissociation in S2B Fig, AN or S1 cultures were incubated for 2 hours 
at 37°C in MB, artificial sea water (ASW) (Instant Ocean, 36 g/L) with different salt concentra- 
tions (18 g/L for 0.5X and 72 g/L for 2X) to assess any effect of electrostatic forces, phosphate- 
buffered saline (PBS) 1X (Sigma-Aldrich) with either Proteinase K at 200 g/mL final concen- 
tration (New England Biolabs, Ipswich, MA, USA) to assess for any protein-dependent effect 
or sonication using a Branson 450 Digital Sonifier (3 pulses of 15 seconds, 10% amplitude). 
Only sonication resulted in dissociation of the clumps and a reduction in sedimentation rates. 
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Maximum likelihood estimation of sedimentation velocity and cellular 
density based on OD600 sedimentation rate assay 


To briefly summarize, we related our OD600 A.U. sedimentation rate measurements and 
radius measurements (S1 Data) to previously published datasets [37] and [91] (S2 and S3 
Data) to gain an estimate of our sedimentation rate measurements in metric units, which 
allowed maximum likelihood estimation of cellular densities based on these estimates 
(S1 Data). 

We started by estimating the average radius of cells from perimeter measurements as 


R = Perimeter/(2 - pi) 
We estimated sedimentation velocity in our measurements by assuming OD600 (OD) 
changes at a constant proportional rate of change with respect to time t = 0 and t = 1 by 


log(OD_1/OD_0) 
t1—t0 





dODdt = 


This yielded a rate of change in arbitrary distance units per second, which we next sought 
to relate to metric distance units. For this, we turned to the datasets from [40], compiled along 
with data from other studies by [37] in his Appendix Table 1 (our S2 Data). In this table, 39 
observations included joint values of salt percentages (or seawater density) for measurement, 
temperature, cell diameters, and sedimentation rate measurements for various phytoplankton 
isolates. 

For each of these observations, we calculated seawater density in experimental assays as SC 
+ salinity where salinity constant SC = 35.16504/35 and the salinity is the salt content in g-I~7 
[91]. We estimated the specific gravity or density of media used in each sedimentation rate 
measurement based upon the dataset published by [91] (S3 Data), using a second-order poly- 
nomial function of their density measurements as a function of salinity and temperature using 
R’s Im() function. The error of this estimate is extremely low (<4e-3 kg-m~*) and was not 
propagated downstream. 

We next calculated the density p, of cells across observations based on sedimentation veloc- 
ity Vin m-s, cell radius R in m, and media density pyby rearrangement of the terminal velocity 
equation 


9-u-V 
Pp = Pf + 2-g-R 
where dynamic viscosity u = 0.00109 Pa-s and gravitational acceleration g = 9.780 m-s*. This 
yielded a median phytoplankton excess density (p,- 1,000) of 139 kg-m”° with a range of 30 to 

1,300 kg-m~*. Some excess density estimates exceeding protein (220 kg-m™) and cellulose (500 
kg-m~*) have previously been suggested to arise by calcide in diatoms [92]. We concluded that 
this method of estimating cell density yielded similar values to those published by, for example, 
[92] and [40]. For downstream analyses, we excluded obvious outliers, including measure- 
ments of samples with densities exceeding that of cellulose, and measurements of D. rex, a 
large diatom, which, in the Smayda dataset, exhibited extraordinarily low sedimentation rates 
given their size. 

To estimate the sedimentation velocities of our dataset, we assumed that our data (S1 Data) 
would fall within the typical measurement in the Smayda dataset (S2 Data). For the outlier- 
excluded subset of the Smayda dataset, we calculated an expected sedimentation velocity for 
what we would measure in our experimental setup based on the specific gravity of the seawater 
formulation we used in our measurements (37.4 g/l or a solvent density of 1,028.9 kg-m~*). We 
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then used least-squares minimization to estimate 2 parameters: a scalar S of the AU velocity 
measurements, dODdt, which could best match our dataset with these expected velocities, and 
an average density parameter Pp par» Which could predict these values and the Smayda dataset’s 
values based on Stokes’ law. For this, we minimized the loss function 


N 


. . . 2 2 
D (vo exp (SH) -* PEE) 





i 


where Vo is either the seawater density-corrected velocity from Smayda or the dODdt parame- 
ter we calculated above, and H is the one-hot binary scalar in [0,1] corresponding to whether 
the 7th data point in our N observations was from Smayda’s or our dataset, respectively. We 
used R’s optim() function with the “L-BFGS-B” method with initial parameters values of 

S = 0.03 and pp nat = 100 and repeated this fit for 500 10-fold (10% out-of-bag) bootstrap sam- 
ples of individual observations across our full dataset and Smayda’s to gain an estimate of the 
error on the parameters. 

S and Pphat fits were positively correlated across testing/training folds (Pearson’s p = 0.97); 
however, mean values were fairly limited, with S = 0.028 + 0.002 and pp nar = 78 + 6.4 kgm? 
(with error equal to the standard deviation of estimates across 500 folds). These narrow esti- 
mates indicated that the fits were reasonably well defined by the underlying dataset. The 
means and standard deviations across predictions of out-of-bag samples are what we have 
reported as means and standard error in S1 Data and propagated along with interreplicate and 
batch error reported in the figures and text. 


Percentage of clumps measurements 


To measure percentage of clumps across all mutants and transfers in S2C Fig, we first soni- 
cated saturated S. arctica cultures using a Branson 450 Digital Sonifier (3 pulses of 15 seconds, 
10% amplitude). Cells were then diluted in fresh MB (1:250 dilution), and cell concentration 
was then measured using a hemocytometer. Approximately 20 cells were then transferred per 
well in a 96-well plate. Cells were then monitored every 24 hours using an inverted optical 
microscope, and the percentage of clumps observed after cellularization was measured manu- 
ally using a tally counter. This experiment was performed 3 independent times, and error bars 
are standard deviations. 


Head-to-head competition 


To perform the fitness experiment (head-to-head competition), we used saturated and soni- 
cated cultures of the ancestor (AN), S01 and S03 evolved isolates. Following the measurement 
of the cell concentration using a neubauer chamber, we independently mixed both S01 and S03 
with the AN at an approximately 1:1 ratio and subjected them to either the same selective 
regime as the evolution experiment (see above) or a control culture that was simply maintained 
in exponential growth without sedimentation selection. The experiment was maintained for 72 
hours with passages every 24 hours. To quantify the outcome of the experiment, every 24 hours 
for 72 hours (7.5 to 12.5 generations), after dilution of the competition culture, we counted the 
number of clumpy and nonclumpy cells at each passage. For this, a fraction of the cultures was 
sonicated and highly diluted to obtain fewer than 10 cells per well of a 96-well plate. These cells 
were allowed to undergo a full life cycle in order to count by microscopy the proportion of cells, 
which bore the clumpy phenotype. The high penetrance of the clumpy phenotype in these 
clones therefore allowed not only a quantification of the clumpy phenotype itself in competi- 
tion, but also the measurement of overall clonal fitness of these derived isolates. 





PLOS Biology | https://doi.org/10.1371/journal.pbio.3001551 March 29, 2022 17/31 


PLOS BIOLOGY 


Regulation of sedimentation rate shapes the evolution of multicellularity in a close relative of animals 





Each data point for the fitness assay was a count of cells with an associated binomial error. 
In Figs 2E and S2D, this error is represented as 95% confidence intervals calculated using the 
binom.logit() function from the R package “binom” (Sundar Dorai-Raj (2014). binom: Bino- 
mial Confidence Intervals for Several Parameterizations. R package version 1.1-1. https:// 
CRAN.R-project.org/package=binom). 

The fitness measurements began from highly synchronized cultures, and, therefore, report- 
ing a per-generation fitness value was clouded by a sudden jump in cellular counts that 
occurred at 48 hours. Therefore, for simplicity and given the clear time-dependent changes in 
clumpy phenotype frequencies, we report fitness in terms of log2 fractional change per unit 
time (“log2 selection coefficient” in days). To obtain these numbers, we first used R’s rbinom() 
function to resample the dataset 500 times, with each resampling using each sample’s original 
measurement of clumpy cell count, total cell count, and proportion clumpy. For each resam- 
pling, after calculating the new value for proportion clumpy, we logit2 transformed the data 
and rescaled to the initial measurement with the following equation: 


Pix Pwo 
lo 2( ) lo 2( ) 
. 1—p,, . 1— Py 


where p is the proportion of the population that was clumpy, the tx underscore is each time 
point, and the t0 underscore is the time point of the initial proportion measurement in units 
days. The slope of the least-squares fit for this quantity with respect to time was obtained by 
R’s lm() function, yielding the per day change for each resampled biological replicate (the 
“log2 selection coefficient”). The mean and standard deviation of this coefficient across all 
resampled replicates were taken as the fitness point estimate and standard error, respectively, 
and the 95% confidence intervals for Fig 2D is that standard error times 1.96. 





Heritability measurements 


To obtain a coarse understanding of how heritable our phenotypes were, we used a statistical 
genetic approach to quantify the heritability of traits, or what proportion of the variance in 
phenotype is due to an individual’s inherited state. Briefly, the total Variance of a measured 
phenotype V_P(Total) can be variance partitioned between Environment, Genetics (or other 
inherited factors), and technical error or stochastic events (epsilon): 


Vor) = Vow + Vive) + Vere) 


The different partitioning terms in V_P(T) can be expressed as a proportion of V_P(T): 


_ VE) + V,(G) + V,(6) 
Vor) 





And, therefore, we can refer to the proportion of variance due to genetics as pV_P(G). 

In the lab, we eliminated environmental variance with tightly controlled experimental con- 
ditions such as temperature and the number of hours a measurement was taken after cell cycle 
synchronization. This allowed us to define heritability therefore as the environment-controlled 
genotypic variance: 


H = pVoein = 1— PV 
For this, we determined H for DNA content, perimeter length, sedimentation rate, and fit- 


ness phenotypes by calculating the variance explained by mean phenotype values within dis- 
tinct genotype/environment combinations (S1A Fig). The results show that H exceeds 95% 
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across phenotypes, and across the entire dataset, H exceeded 99% of the total phenotypic vari- 
ance (ANOVA F = 1,118 on 252 and 735 DF, p = 0). This means that for a typical individual 
genotype in a given environment, we could predict its average phenotypic measurement with 
>97% accuracy. 


Microscopy 


Microscopy of live and fixed cells was performed using a Zeiss Axio Observer Z.1 Epifluores- 
cence inverted microscope equipped with Colibri LED illumination system and an Axiocam 
503 mono camera. An EC Plan-Neofluar 40x/0.75 air objective was used for images of fixed 
cells, and an N-ACHROPLAN 20x/0.45na Ph2 air objective was used for all live imaging, 
unless indicated otherwise. 


Cell fixation and staining 


Throughout this study, saturated Sphaeroforma cultures were mildly sonicated prior to dilut- 
ing them 250X in fresh marine broth to initiate a synchronized culture. To assess for any tem- 
perature dependency, cultures were grown at both 17°C and 12°C, and measurements were 
conducted every 12 hours for a duration of 72 hours. For every time point, cells were fixed 
using 4% formaldehyde and 250 mM sorbitol for 30 minutes before being washed twice with 
PBS. For nuclei staining, cells were centrifuged at 1,000 rpm for 3 minutes after fixation and 
washed again 3 times with PBS before adding DAPI at a final concentration of 5 ug/mL to 5 ul 
of concentrated sample. DAPI-stained samples were imaged to measure DNA content and 
coenocyte size. It is important to note that results obtained from fast-settling mutants prior to 
cell release correspond to measurements of unicellular coenocytes (24 hours for 17°C and 36 
hours for 12°C), whereas results collected after cell release correspond to measurements of 
multicelled clumps (48 hours for 17°C and 72 hours for 12°C). For cell wall staining, cells were 
incubated with Calcofluor-white (Sigma-Aldrich) at a final concentration of 5 ug/ml from a 
200X stock solution prior to fixation. Cells were then fixed as previously mentioned and con- 
centrated before being disposed between slide and coverslip. 


Live cell imaging 

For live cell imaging, saturated cultures were diluted 250x in fresh MB medium inside a - 
Slide 4- or 8-well slide (Ibidi) at time 0. To ensure oxygenation during the whole period of the 
experiment, the cover was removed. To maintain the temperature at 17 or 12°C, we used a 
P-Lab Tek (Pecon GmbH) Heating/Cooling system connected to a Lauda Ecoline E100 circu- 
lating water bath. To reduce light toxicity, we used a 495-nm Long Pass Filter (FGL495M- 
ThorLabs). For plasma membrane live staining (Fig 3B, S3 Video), FM4-64 (Invitrogen) at a 
final concentration of 10 uM from a 100x DMSO-diluted stock solution was added at time 0 
unless indicated otherwise in figure legends. For cytoplasmic staining of cells in S3B Fig, cells 
were either stained with CellTrace CFSE Cell Proliferation Kit (Thermofisher) or CellTrace 
Calcein Red-Orange (Thermofisher). 


Image analysis 


Image analysis was done using Image] software (version 1.52) [93]. For nuclear content distri- 
bution across Sphaeroforma sp.’s life cycle, fixed and DAPI-stained coenocytes were imaged, 
and the number of nuclei per coenocyte was counted using the Object] plugin in image]. To 
compute nuclear duplication times, log2 of geometric mean of DNA content was calculated as: 
log,(geommean) = 7 f, « log, (x;) where fi is the fraction of cells and xi the DNA content 
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(number of nuclei per cell) (ploidy) of each i-th DNA content bin. Nuclear doubling times 
were computed as linear regression of log2 of geometric mean of DNA content versus time. 
Note that for S. tapetis, nuclear doubling times could not be computed due to the asynchrony 
in growth. For measurements of cell volume in live and fixed cells, we used the oval selection 
tool to draw the contour of each cell and measured cell perimeter. As cells are spherical, we 
computed cell volume as V = 4/3zr° where r is the cell radius. For measurements of clumps 
perimeter, we transformed the images into binaries to ensure later segmentation. We then 
used the particle analysis function in Image] with a circularity parameter set to 0.15 to 1 to 
measure cell perimeter. For nuclear number-to-volume ratios, the number of nuclei was 
divided by the coenocyte volume measured as previously described for fixed cells. All figures 
were assembled with Illustrator CC 2020 (Adobe). Several figures were generated using ggplot2 
in R version 4.0.5 [94]. 


Sphaeroforma arctica genome sequencing and assembly 


Genomic DNA was extracted for the ancestral strain (AN) and a single clonal isolate from 
each evolved population (S1-S10) using QlAamp DNA Blood Midi Kit (Qiagen) following the 
manufacturer’s recommendations from 50 mL culture incubated at 17°C for 5 days in 75 cm? 
flasks. The Qubit (Invitrogen) quantification ranged between 3 and 13 ug of genomic DNA in 
total. All of the subsequent steps were performed by the CRG Genomics Unit (Barcelona): 
Sequencing libraries were prepared from the pure high molecular weight DNA using TruSeq 
DNA HT Library Preparation kit (Illumina HiSeq Sequencing v4 Chemistry). A paired-end 
library with a target insert size of approximately 500 bp was sequenced on an Illumina 
HiSeq2500 platform in paired-end mode, with read lengths of 125 bp. The resulting paired raw 
read files were demultiplexed by the sequencing facility and data stored in 2 separate, gzip- 
compressed FASTQ files of equal sizes. Genome sequencing data have been deposited in NCBI 
SRA under the BioProject accession PRJNA693121. 


Bioinformatic analyses of the genomes 


Data processing. On average, each paired-end sequencing library contained approxi- 
mately 58.8 million reads of 125 bp sequence lengths (S2 Table), equaling approximately 7.35 
billion base pairs (Gbp). From these data, we carefully removed adapter sequences and reads 
shorter than 50 bp from the raw read data using trimmomatic v0.36 [95], yielding an average 
of approximately 4.36 Gbp of filtered sequence data per genome (S2 Table). The quality of 
both raw and trimmed sequencing data was assessed in FastQC v0.11.7 [96]. In more detail, 
based on the FastQC output for raw reads, we initiated trimming by calling the following 
parameters: 

ILLUMINACLIP:Nexterat+TruSeq3-PE2.fa:3:25:10 CROP:110 LEADING:30 TRAILING:25 
SLIDINGWINDOW:4:28 MINLEN:50. 

This translates into the following trimming steps: 


« cut adapters and other Ilumina-specific sequences using a combined file of default adapters 
(“Nextera.fa” AND “TruSeq3-PE2.fa”) to catch as many spurious contaminations during 
library prep as possible, with seed mismatches = 3; palindrome clip threshold = 25; simple 
clip threshold = 10; 


« end-clipping of the final 15 bp of all reads, due to evidence for elevated adapter content; 
¢ quality-clip all bases on leading ends as long as bases were of lower quality than Q < 30; 


* removing all bases on trailing ends as long as bases were of lower quality than Q < 25; and, 
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¢ finally, conducting a sliding window approach, where the reads were trimmed once the aver- 
age quality within a window of 4 consecutive bases falls below a threshold of Q <28. 


The trimming output consisted of 4 FASTQ files for each genome, of which 2 files con- 
tained intact paired-end reads, and another 2 files containing all unpaired reads for each end 
separately. 

Repeat-masking of reference genome. We relied on the latest (i.e., fourth) assembly ver- 
sion of the S. arctica reference genome (Sarc4; [32]) for variant detection and annotation. Ini- 
tial runs of the read-mapping steps revealed a high proportion of tightly clustered variants 
(approximately 20.4%) concentrated in certain regions of the genome (S3 Table). We investi- 
gated this phenomenon and determined that the issue was caused by repetitive stretches of 
sequence and thus decided to mask all of the potentially problematic repeat regions. For this, 
repetitive regions were screened for and properly annotated in RepeatMasker v4.1.1 [97] rely- 
ing on the 20181026 release of GIRI RepBase database for annotations [98] and applying the 
slow high-sensitivity search mode with the following parameters: -pa 10 -s -gff -excln -species 
Opisthokonta. 

Variant prediction. We performed read alignment, variant calling, variant filtering in 
CLC Genomic Workbench v20.0.4 (Qiagen), and our analytical pipeline was structured into 
the following steps: 


¢ Read alignment. Paired-end reads were merged (assuming insert lengths between 400 and 
600 bp), and both paired and unpaired reads were aligned against the reference genome 
(with settings: Length fraction > 0.9; Similarity fraction > 0.9; Match score = 1; Mismatch 
cost = 3; Insertion/Deletion Open Cost = 5; Insertion/Deletion extend cost = 3; Global 
Alignment = no). Finally, reads were deduplicated (maximum representation of minority 
sequence = 0.2). 


2 


¢ Variant calling. We called variants using CLC’s “Fixed Ploidy Variant Detection” assuming a 
haploid genome (Ploidy = 1) for S. arctica. We further required a variant probability of at 
least 90%, ignoring positions in excess of 90x coverage, broken read pairs and nonspecific 
matches. All variants needed to be covered by at least 10 variant-bearing reads and a mini- 
mum consensus of 80% (i.e., at least 10/12 variant supporting reads). We also applied the fol- 
lowing read quality filters: neighborhood radius = 10; minimum central quality = 20; 
minimum neighborhood quality = 20), and read direction filters (direction frequency = 0.05; 
relative read direction filter = yes (significance = 0.05); read position filter = yes 
(significance = 0.05). 


Variant filtering. (1) To consider only variants that have emerged through the course of evo- 
lution, we automatically removed mutations present in the ancestor (AN). (2) We then went 
on to manually curate all mutation predictions. We did this by aligning the mapping tracks 
(profiles) of all resequenced genomes (AN, S1-S10) and screening all initial candidate muta- 
tions. The overwhelming majority of variant calls were visually shared across all evolved 
clones but not called universally due to extremely low statistical support, or low local cover- 
age in some of the genomes. The final dataset hence only contained 25 predicted variants 
(one shared among 2 lineages, hence, representing 26 independent mutations) and exported 
these for each of the 10 evolved single clones as Variant Call File (VCF) format. 


¢ Variant annotation. We annotated filtered variants from converted VCF files with breseq 
v0.33.2 [99] for convenience. More specifically, we converted VCF files into breseq’s Genome 
Diff (GD) file format using the command “gdtools VCF2GD”. We then annotated all 
genomes in GD format jointly by running the command “gdtools ANNOTATE” and 








PLOS Biology | https://doi.org/10.1371/journal.pbio.3001551 March 29, 2022 21/31 


PLOS BIOLOGY 


Regulation of sedimentation rate shapes the evolution of multicellularity in a close relative of animals 





specifying original S. arctica genome assembly (Sarc4) in GenBankFormat (GBK) as refer- 
ence. Mutations were tabulated and sorted. Finally, we counted and categorized both 
observed mutations and nonsynonymous and synonymous sites at risk across the reference 
genome using the command “gdtools COUNT -b” for statistics and the calculation of the 
compound dN/dS ratio. 


Supporting information 


$1 Data. Our data used for estimating sedimentation rate in metric units and density in 
mass per volume. “V_au” is sedimentation rate in AU OD600 units per second. “R_meters” is 
the mean radius of cells for the genotype and environment. “Species,” genotype, temp, and 
hours_growth are information for aggregating the data based on genotype and environment. 
“V_mu” is the maximum likelihood estimate of cellular sedimentation rate—the mean of all 
out-of-bag samples used in fitting—in meters per second. “V_se” is the standard deviation of 
cellular sedimentation rate across all out-of-bag samples in the fitting. “pp_mu” the maximum 
likelihood estimate of density—the mean of all out-of-bag samples used in fitting—in kg/m’. 
“pp_se” the maximum likelihood estimate of density—the standard deviation of all out-of-bag 
samples used in fitting—in kg/m’, 

(XLSx) 


$2 Data. Smayda dataset used for calibrating our dataset. Adapted from [37] Appendix 
Table 1. “classification” species are from his table annotation. “salt_percent” is the percent salt 
reported in the table. “V_meters_per_second” is the velocity reported, converted from meters 
per day to meters per second. “Pf” is the density of the media based on the reported tempera- 
ture and salt concentration, estimated by the data from Millero and Huang. “R_meters” is the 
mean radius in meters. When a range was given, this is the average of that range. “Pp_” is the 
density of the sample in kg/m’. “Pp” is the excess density kg/m*. “V_exp” is the expected veloc- 
ity in the seawater used in our experiments (37.4 g/l). 


(XLSX) 


$3 Data. The Millero and Huang data [92] used for estimating seawater density based on 
salinity and temperature. “temp_C” temperature in degrees Celsius. “salinity” salt in g/l. 
“density_kgm3” excess density in kg/m’. “density” density in kg/m’. 

(XLSX) 


S4 Data. Fitness experiment data used for heritability measurements. “batch” is an arbitrary 
alphabetical letter referring to the date the replicate was performed (either D or E). “rep” is the 
replicate within batch/day (rep1 or rep2). “geno” is the genotype information (which genotype 
was being competed against the “AN” ancestor, or “alone” if alone in competition). “time” is 
the time the sample was collected in units’ hours. “sel_cond” is the condition of the selection— 
either the initial population (“init”), no selection for sedimentation rate (“NO_SEL”), selection 
for sedimentation rate (“YES_SEL”). “NOT_CLUMPY” is the count of microcolonies with not 
clumpy phenotype. “YES_CLUMPY” is the count of microcolonies with clumpy phenotype. 
“cell_density” is the measured cell density at time of dilution after period of growth. “N” is the 
total population growth calculated from cell_density per fold dilutions (1/250 per day). “gen” 
is the number of generations of growth based on cellular counts. Note for both “N” and “gen” 
that due to the synchronized coenocytic growth the cultures, they underwent a large increase 
in cell count after 48 hours. “date” is the date the replicate was performed (April or May 2019). 
“temp” is the temperature of the competition. The data underlying this figure may be found in 
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S4 Data; see Methods. 
(XLSX) 


$1 Fig. Variance measurements and heritability of the characterized phenotypes across 
this study. (A) Measurements of technical error or heritability of traits show that the pheno- 
types measured for the evolved isolates were highly heritable. Arithmetic mean values for 4 key 
phenotypes were calculated for each genotype + environment combination (horizontal X axis) 
and plotted against individual measured values (vertical Y axis) across replicates performed on 
separate days (“batches”—colors). Environments included different temperatures and different 
timings after cell cycle synchronization. Horizontal facet strips indicate the phenotype, its unit 
of measure, and the heritability of the trait, i-e., the variance explained by the mean genotype 

+ environment value. Gray lines represent where X = Y. 

(TIF) 


$2 Fig. Clump size of evolved mutants correlates with increased number of cells per 
clump. (A) Approximation of the constant population size derived from nuclear doublings 
over the course of 48 hours. Nuclei numbers were estimated from total cell number and aver- 
age number of nuclei per cells (after cell sorting). At 24 hours, a constant fraction of cells has 
been transferred to fresh medium, which induces a bottleneck. The data underlying this figure 
may be found in S1 Table. (B) Sedimentation rates of S1 clumps after incubation for 2 hours in 
different media or after sonication. MB = marine broth, ASW = artificial sea water (Salt con- 
centration = 36.4 g/L), 2X (Salts concentration = 72.8 g/L), 0.5X (Salts concentration = 18.2 g/ 
L), PBS = phosphate-buffered saline 1X, Proteinase K 200 ug/mL and sonication (3 pulses of 
15 seconds, 10% amplitude). Representative images on the right show the dissociation of $1 
clumps after sonication. Bar, 50 um. (C) Percentage of clumps formed after 72 hours of all 
evolved mutants per number of transfers (n > 200 coenocytes at cell release per strain). Exam- 
ple images on the left of detached or clumpy cells after cell release. Bar, 50 um. (D) Raw data 
for the frequency with respect to time of each independent replicate of strains ancestor (AN) 
and clumpy evolved isolates S01 and S03 in head-to-head competitions with the ancestral 
strain (2 days, 2 replicates per day). Each point represents a measurement of a given replicate 
at the given time point under the 2 selection regimes (colors—Green = in presence of sedimen- 
tation rate selection; Black = in absence of selection), and error bars are 95% confidence inter- 
vals. Lines are least-squares regressions of each replicate across time points to guide the eye. 
The data underlying this graph may be found in S4 data; see Methods. (E) Number of cells per 
clump in all evolved mutants illustrates the linear correlation between sedimentation rate and 
number of cells per clump. (” = 13 clumps for AN and >50 clumps for all evolved mutants of 
at least 5 attached cells together). (F) Linear correlation between clump size of the fast-settling 
mutants (S1, S4, and $9) and number of cells per clump compared to the ancestral strain (AN). 
The data underlying this figure may be found at https://doi.org/10.6084/m9.figshare. 18319232. 
vl. 

(TIF) 


$3 Fig. Fast-settling mutants cannot form clumps by aggregation and show discreet differ- 
ences in life stages duration. (A) Cell wall staining show presence of a separating cell wall 
between individual cells in the clumps of all fast-settling mutants. Bar, 50 um. (B) Experimen- 
tal design and measurements of S1 clump formation by aggregation. AN and S1 cells are sepa- 
rated by sonication and stained with different cellular dyes before being mixed together for a 
complete life cycle of 72 hours. Clumps of each separate or mixed colors are then counted by 
microscopy. Circles represent the number of clumps observed in each condition (more than 5 
cells attached together). We show that when S1 is mixed with either AN or itself (1:1 ratio), it 
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mostly forms monocolored clumps. Representative images of monocolored clumps. Bar, 

50 um. (C) Mean coenocyte perimeter over time (10 cell traces per strain) at 12°C, aligned to 
time 0, reveals discreet differences in coenocyte perimeter and life cell stages among fast-set- 
tling mutants. Orange squares represents the flip time point in each trace. (D) Duration of 
growth, cellularization, and post-flip represented as box-plots at 12°C (n > 28 coenocytes 
each). The data underlying this figure may be found at https://doi.org/10.6084/m9.figshare. 
183192321. 

(TIF) 


$4 Fig. Increased sedimentation rates in fast-settling mutants is independent of tempera- 
ture. (A) Sedimentation rates of S. arctica AN and evolved mutants during the life cycle at 
12°C. Every trace represents an independent experiment. (B) Average perimeter measured 
from fixed cells every 12 hours over a complete life cycle of 72 hours at 12°C shows that fast- 
settling mutant increase their size upon cell release. Every trace represents an independent 
experiment (n > 130 measurements per time point for each independent experiment). (C) 
Average perimeter of fast-settling cells and clumps at 36 and 72 hours, respectively, show that 
S4 and S9 single cells and clumps have a smaller size when compared to S1 at 12°C. Every 
square represents an independent experiment, and the white circle represents the median. (D) 
Excess cellular density of fast-settling mutants (before cellularization) and clumps (after cellu- 
larization) at 36 and 72 hours, respectively. Every square represents an independent experi- 
ment, and the white circle represents the median. (E) Distributions of nuclear content of S. 
arctica AN and fast-settling mutants during the life cycle at 17°C measured by microscopy of 
DAPI-fixed cells (n > 400 coenocytes per time point). (F) Distributions of nuclear content of 
S. arctica AN and fast-settling mutants during the life cycle at 12°C measured by microscopy 
of DAPI-fixed cells (n > 420 coenocytes per time point for each independent experiment). (G) 
Quantification of mean DNA content per time point for fast-settling mutants grown in marine 
broth at 12°C. Every trace represents an independent experiment (n > 420 coenocytes per 
time point for each independent experiment). (H) Nuclear doubling time, calculated by linear 
regression of mean nuclear content at time points from 0 to 24 hours at 12°C. Every square 
represents an independent experiment, and the white circle represents the median (n > 420 
coenocytes per time point for each independent experiment). (I) Boxplots of cell volume mea- 
surements of DAPI-stained fixed cells at 12°C. For 1-, 4-, 16-, and 64-nuclei cells. Cells with 1 
nucleus represent newborn cells at the end of the experiment (n > 50 coenocytes per DNA 
content). (J) Boxplots of nuclear number-to-volume ratio of DAPI-stained cells at 12°C. Every 
square represents an independent experiment, and the white circle represents the median 

(n > 300 coenocytes per strain). (K) Temporal transcript abundance of genes mutated in inter- 
mediate and slow-settling phenotypes across the native life cycle of S. arctica. The data under- 
lying these graphs may be found in S5 Table. The data underlying this figure may be found at 
https://doi.org/10.6084/m9.figshare.18319232.v1. 

(TIF) 


$5 Fig. Sphaeroforma sister species are distinct in growth, sedimentation rates, and cell 
volume. (A) Distributions of nuclear content of Sphaeroforma species during the life cycle at 
17°C of DAPI-fixed cells measured by microscopy. Note that S. tapetis is asynchronous com- 
pared to all other species (n > 300 coenocytes per time point for each independent experi- 
ment). (B) Average perimeter measured from fixed cells every 12 hours over a complete life 
cycle of 72 hours at 17°C of Sphaeroforma sp. Every trace represents an independent experi- 
ment (n > 140 measurements per time point for each independent experiment). (C) Average 
perimeter of Sphaeroforma sp. cells and clumps at 0, 24, and 48 hours, respectively, show that 
S. gastrica, S. nootakensis, S. napiecek, and S. sirkka are smaller in size at 17°C. Every square 
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represents an independent experiment, and the white circle represents the median (n > 140 
measurements per time point for each independent experiment). (D) Coenocyte volume at 
flip, measured from time-lapse movies, show that all Sphaeroforma species apart from S. arc- 
tica have significantly smaller coenocyte volume at flip (n > 50 coenocytes per strain). (E) Box- 
plots of cell volume measurements of DAPI-stained fixed Sphaeroforma sp. coenocytes at 
17°C. For 1-, 4-, 16-, and 64-nuclei cells. Cells with 1 nucleus represent newborn cells at the 
end of the experiment. (n > 50 coenocytes per DNA content). The data underlying this figure 
may be found at https://doi.org/10.6084/m9.figshare.18319232.v1. 

(TIF) 


S1 Video. Video of Sphaeroforma arctica AN and fast-settling (S1, $4, and S9) cultures sed- 
imenting. Time interval between frames is 0.5 second. The movie is played at 7 frames per sec- 
ond (fps). We can observe rapid cell sedimentation in S1, $4, and S9 when compared to AN. 
Note that S1 clumps are bigger and sediment faster than the 2 other mutants. The movie was 
acquired for cultures pregrown for 72 hours at 12°C and obtained with a mobile phone (Sam- 
sung A20). 

(MP4) 


$2 Video. Time-lapse video of synchronized cells of S. arctica AN and fast-settling mutants 
(S1, $4, and $9). Time interval between frames is 30 minutes. The movie is played at 7 fps. 
Four distinct cells can be seen undergoing a full life cycle at 12°C with the release of detached 
newborn cells for AN or clumps for the mutants. Bar, 50 tm. 

(MP4) 


S3 Video. Time-lapse video of cells of S. arctica AN and fast-settling mutants (S1, $4, and 
S9) stained with the membrane dye (Fm4-64). Time interval between frames is 15 minutes. 
The movie is played at 7 fps. Clumps can be seen being formed after plasma membrane invagi- 
nations followed by cell release. Obtained with an epifluorescence microscope. Bar, 50 um. 
(MP4) 


$4 Video. Time-lapse video of 6 different Sphaeroforma species undergoing a complete life 
cycle. Time interval between frames is 30 minutes. The movie is played at 7 fps. Note the asyn- 
chrony of S. tapetis, the capacity to clump of S. gastrica and S. nootakensis, and the small cell 
size of all Sphaeroforma species compared to S. arctica. Bar, 50 um. 

(MP4) 


$1 Table. Inference of doubling times and effective population size [Ne] from direct esti- 
mates of cell numbers and nuclei per cells over the time course of 2 transfers. 
(PDF) 


$2 Table. Genome resequencing summary: Data processing and mapping. 
(PDF) 


$3 Table. Summary and classification of annotated repeat regions in the Sphaeroforma arc- 
tica reference genome. 
(PDF) 


$4 Table. Variant calling. 
(PDF) 


$5 Table. Transcript abundance of affected genes during S. arctica native life cycle. 
(PDF) 
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S6 Table. Annotation and gene prediction of mutational targets. 
(PDF) 


Acknowledgments 


We thank Pierre Génczy, Macarena Toll Riera, Gautam Dey, and Andrej Ondracka for discus- 
sion and comments on the manuscript; Hiroshi Suga and Xavi Grau-Bové for advice on 
genome analysis; Jon Brate for cultures of the different Sphaeroforma sister species; and Merit- 
xell Anto for technical support. We also thank Andrej Ondracka for sharing his unpublished 
results on Sphaeroforma sister species and Thibaut Brunet for sharing his unpublished obser- 
vations on the sedimentation phenotype of Choanoflagellate colonies. We also acknowledge 
the CRG Genomics Unit for mRNA library preparation and Illumina sequencing. 


Author Contributions 

Conceptualization: Omaya Dudin. 

Data curation: Omaya Dudin, Sébastien Wielgoss, Aaron M. New. 
Formal analysis: Omaya Dudin, Sébastien Wielgoss, Aaron M. New. 
Funding acquisition: Omaya Dudin, Ifaki Ruiz-Trillo. 
Investigation: Omaya Dudin, Sébastien Wielgoss, Aaron M. New. 
Methodology: Omaya Dudin, Sébastien Wielgoss, Aaron M. New. 
Project administration: Omaya Dudin. 

Supervision: Ifaki Ruiz-Trillo. 

Visualization: Omaya Dudin. 

Writing - original draft: Omaya Dudin, Sébastien Wielgoss, Aaron M. New. 


Writing - review & editing: Omaya Dudin, Sébastien Wielgoss, Aaron M. New, Ifaki Ruiz- 
Trillo. 


References 


1. Parfrey LW, Lahr DJG. Multicellularity arose several times in the evolution of eukaryotes. BioEssays. 
2013. https://doi.org/10.1002/bies.201200143 PMID: 23315654 


2. Sebé-Pedrdés A, Degnan BM, Ruiz-Trillo |. The origin of Metazoa: A unicellular perspective. Nat Rev 
Genet. 2017. https://doi.org/10.1038/nrg.2017.21 PMID: 28479598 


3. Grosberg RK, Strathmann RR. The evolution of multicellularity: A minor major transition? Annu Rev 
Ecol Evol Syst. 2007. https://doi.org/10.1146/annurev.ecolsys.38.091206.095634 PMID: 25685107 


4. Knoll AH. The Multiple Origins of Complex Multicellularity. Annu Rev Earth Planet Sci. [Internet]. 2011 
May 30 [cited 2021 Feb 23]; 39(1):217-39. Available from: http://www.annualreviews.org/doi/10.1146/ 
annurev.earth.031208.100209. 


5. Niklas KJ, Newman SA. The origins of multicellular organisms. Evol Dev [Internet]. 2013 Jan [cited 
2021 Feb 23]; 15(1):41—-52. Available from: https://pubmed.ncbi.nim.nih.gov/23331916/. https://doi.org/ 
10.1111/ede.12013 PMID: 23331916 


6. Rokas A. The molecular origins of multicellular transitions. Curr Opin Genet Dev. 2008; 18:472-8. 
https://doi.org/10.1016/j.gde.2008.09.004 PMID: 18926910 


7. Ruiz-Trillo |, Burger G, Holland PWH, King N, Lang BF, Roger Au, et al. The origins of multicellularity: a 
multi-taxon genome initiative. Trends Genet [Internet]. 2007 Mar [cited 2021 Feb 23]; 23(3):113-8. 
Available from: https://pubmed.ncbi.nlm.nih.gov/17275133/. https://doi.org/10.1016/j.tig.2007.01.005 
PMID: 17275133 


8. Leigh EG, Smith JM, Szathmary E. The Major Transitions of Evolution. Evolution (N Y). 1995. 





PLOS Biology | https://doi.org/10.1371/journal.pbio.3001551 March 29, 2022 26/31 


PLOS BIOLOGY 


Regulation of sedimentation rate shapes the evolution of multicellularity in a close relative of animals 





17. 


18. 


19. 


20. 


21. 


22. 


23. 


24. 


25. 


26. 


27. 


28. 


29. 


30. 


31. 


Bonner JT. The origins of multicellularity. Integr Biol [Internet]. 1998 Jan 1 [cited 2021 Feb 23]; 1(1):27- 
36. Available from: https://onlinelibrary.wiley.com/doi/10.1002/(SICI)1520-6602(1998)1:1%3C27::AID- 
INBI4%3E3.0.CO;2-6. 


Tarnita CE, Taubes CH, Nowak MA. Evolutionary construction by staying together and coming 
together. J Theor Biol. 2013. https://doi.org/10.1016/j.jtbi.2012.11.022 PMID: 23206384 


Wielgoss S, Wolfensberger R, Sun L, Fiegna F, Velicer GJ. Social genes are selection hotspots in kin 
groups of a soil microbe. Science. 2019. https://doi.org/10.1126/science.aar4416 PMID: 30898932 


Kawecki TJ, Lenski RE, Ebert D, Hollis B, Olivieri 1, Whitlock MC. Experimental evolution. Trends Ecol. 
2012. 


Elena SF, Lenski RE. Evolution experiments with microorganisms: the dynamics and genetic bases of 
adaptation. Nat Rev Genet [Internet]. 2003; 4(6):457-69. Available from: https://doi.org/10.1038/ 
nrg1088 PMID: 12776215 


Ratcliff WC, Fankhauser JD, Rogers DW, Greig D, Travisano M. Origins of multicellular evolvability in 
snowflake yeast. Nat Commun. 2015. https://doi.org/10.1038/ncomms7 102 PMID: 25600558 


Ratcliff WC, Denison RF, Borrello M, Travisano M. Experimental evolution of multicellularity. Proc Natl 
Acad Sci US A. 2012. https://doi.org/10.1073/pnas. 1115323109 PMID: 22307617 


Ratcliff WC, Herron MD, Howell K, Pentz JT, Rosenzweig F, Travisano M. Experimental evolution of an 
alternating uni- and multicellular life cycle in Chlamydomonas reinhardtii. Nat Commun. 2013. https:// 
doi.org/10.1038/ncomms3742 PMID: 24193369 


Koschwanez JH, Foster KR, Murray AW. Improved use of a public good selects for the evolution of 
undifferentiated multicellularity. elife. 2013 Apr 2;2013(2). https://doi.org/10.7554/eLife.00367 PMID: 
23577233 


Herron MD, Borin JM, Boswell JC, Walker J, Chen ICK, Knox CA, et al. De novo origins of multicellular- 
ity in response to predation. Sci Rep. 2019. https://doi.org/10.1038/s41598-019-39558-8 PMID: 
30787483 


Oud B, Guadalupe-Medina V, Nijkamp JF, De Ridder D, Pronk JT, Van Maris AJA, et al. Genome dupli- 
cation and mutations in ACE2 cause multicellular, fast-sedimenting phenotypes in evolved Saccharo- 
myces cerevisiae. Proc Natl Acad Sci U S A [Internet]. 2013 Nov 5 [cited 2021 Feb 23]; 110(45):E4223- 
31. Available from: www.pnas.org/cgi/doi/10.1073/pnas.13059491 10. PMID: 24145419 


Strassmann JE, Gilbert OM, Queller DC. Kin discrimination and cooperation in microbes. Annu Rev 
Microbiol. 2011. https://doi.org/10.1146/annurev.micro.112408.134109 PMID: 21682642 


Fisher RM, Cornwallis CK, West SA. Group formation, relatedness, and the evolution of multicellularity. 
Curr Biol. 2013 Jun 17; 23(12):1120—5. https://doi.org/10.1016/j.cub.2013.05.004 PMID: 23746639 


Queller DC, Ponte E, Bozzaro S, Strassmann JE. Single-gene greenbeard effects in the social amoeba 
Dictyostelium discoideum. Science. 2003. https://doi.org/10.1126/science. 1077742 PMID: 12511650 


Brunet T, King N. The Origin of Animal Multicellularity and Cell Differentiation. Dev Cell. 2017. https:// 
doi.org/10.1016/j.devcel.2017.09.016 PMID: 29065305 


Ruiz-Trillo |, de Mendoza A. Towards understanding the origin of animal development. Development. 
2020. https://doi.org/10.1242/dev. 192575 PMID: 33272929 


Pérez-Posada A, Dudin O, Ocafia-Pallarés E, Ruiz-Trillo |, Ondracka A. Cell cycle transcriptomics of 
Capsaspora provides insights into the evolution of cyclin-CcDK machinery. PLoS Genet. 2020; 16 (3). 
https://doi.org/10.1371/journal.pgen.1008584 PMID: 32176685 


Parra-Acero H, Harcet M, Sanchez-Pons N, Casacuberta E, Brown NH, Dudin O, et al. Integrin-Medi- 
ated Adhesion in the Unicellular Holozoan Capsaspora owczarzaki. Curr Biol. 2020. https://doi.org/10. 
1016/j.cub.2020.08.015 PMID: 32857975 


Brunet T, Larson BT, Linden TA, Vermeij MJA, McDonald K, King N. Light-regulated collective contrac- 
tility in a multicellular choanoflagellate. Science. 2019. 


Levin TC, Greaney AJ, Wetzel L, King N. The Rosetteless gene controls development in the choanofla- 
gellate S. rosetta. elife. 2014. https://doi.org/10.7554/eLife.04070 PMID: 25299189 


Fairclough SR, Dayel Mu, King N. Multicellular development in a choanoflagellate. Curr Biol. 2010. 


Alegado RA, Brown LW, Cao S, Dermenjian RK, Zuzow R, Fairclough SR, et al. A bacterial sulfonolipid 
triggers multicellular development in the closest living relatives of animals. elife. 2012. https://doi.org/ 
10.7554/eLife.00013 PMID: 23066504 


Sebé-Pedrés A, Irimia M, del Campo J, Parra-Acero H, Russ C, Nusbaum C, et al. Regulated aggrega- 
tive multicellularity in a close unicellular relative of metazoa. elife. 2013 Dec;2( 2:e01287). https://doi. 
org/10.7554/eLife.01287 PMID: 24368732 





PLOS Biology | https://doi.org/10.1371/journal.pbio.3001551 March 29, 2022 27/31 


PLOS BIOLOGY 


Regulation of sedimentation rate shapes the evolution of multicellularity in a close relative of animals 





32. 


33. 


34. 


35. 


36. 


37. 


38. 


39. 


40. 


41. 


42. 


43. 


44. 
45. 


46. 


47. 


48. 


49. 


50. 


51. 


52. 


Dudin O, Ondracka A, Grau-Bové X, Haraldsen AAB, Toyoda A, Suga H, et al. A unicellular relative of 
animals generates a layer of polarized cells by actomyosin-dependent cellularization. elife. 2019; 8. 
https://doi.org/10.7554/eLife.49801 PMID: 31647412 


de Mendoza A, Suga H, Permanyer J, Irimia M, Ruiz-Trillo |. Complex transcriptional regulation and 
independent evolution of fungal-like traits in a relative of animals. elife. 2015 Oct; 4. 


Ondracka A, Dudin O, Ruiz-Trillo |. Decoupling of Nuclear Division Cycles and Cell Size during the Coe- 
nocytic Growth of the Ichthyosporean Sphaeroforma arctica. Curr Biol. 2018; 28 (12). https://doi.org/10. 
1016/j.cub.2018.04.074 PMID: 29887314 


Glockling SL, Marshall WL, Gleason FH. Phylogenetic interpretations and ecological potentials of the 
Mesomycetozoea (Ichthyosporea). Vol. 6, Fungal Ecol. Elsevier Ltd; 2013. p. 237-47. 


Mendoza L, Taylor JW, Ajello L. The class Mesomycetozoea: A heterogeneous group of microorgan- 
isms at the animal-fungal boundary [Internet]. Annu Rev Microbiol. 2002 [cited 2021 Feb 23]; 56:315— 
44. Available from: https://pubmed.ncbi.nim.nih.gov/12142489/. https://doi.org/10.1146/annurev.micro. 
56.012302.160950 PMID: 12142489 


SMAYDA JT The suspension and sinking of phytoplankton in the sea. Ocean Mar Biol Ann Rev [Inter- 
net]. 1970 [cited 2021 Apr 13]; 8:353-414. Available from: https://ci.nii.ac.jp/naid/10009425107. 


Beardall J, Allen D, Bragg J, Finkel Z V, Flynn KJ, Quigg A, et al. Allometry and stoichiometry of unicellu- 
lar, colonial and multicellular phytoplankton. New Phytol [Internet]. 2009 Jan 1 [cited 2021 Apr 13]; 181 
(2):295-309. Available from: http://doi.wiley.com/10.1111/j.1469-8137.2008.02660.x. PMID: 19121029 


Finkel Z V., Beardall J, Flynn KJ, Quigg A, Rees TA V, Raven JA. Phytoplankton in a changing world: 
Cell size and elemental stoichiometry [Internet]. J Plankton Res. Oxford Academic; 2010 [cited 2021 
Apr 13]; 32:119-37. Available from: https://academic.oup.com/plankt/article/32/1/119/1492394. 


Eppley RW, Holmes RW, Strickland JDH. Sinking rates of marine phytoplankton measured with a fluo- 
rometer. J Exp Mar Bio Ecol. 1967 Jan 1; 1(2):191-208. 


Smayda TJ. Normal and accelerated sinking of phytoplankton in the sea. Mar Geol. 1971 Sep 1; 11 
(2):105-22. 

Gemmell Bu, Oh G, Buskey Eu, Villareal TA. Dynamic sinking behaviour in marine phytoplankton: 
Rapid changes in buoyancy may aid in nutrient uptake. Proc R Soc B Biol Sci [Internet]. 2016 Oct 12 
[cited 2021 Apr 13]; 283(1840). Available from: https://royalsocietypublishing.org/doi/abs/10.1098/rspb. 
2016.1126. PMID: 27708154 


Friebele ES, Correll DL, Faust MA. Relationship between phytoplankton cell size and the rate of ortho- 
phosphate uptake: in situ observations of an estuarine population. Mar Biol [Internet]. 1978 Mar [cited 
2021 Apr 13]; 45(1):39-52. Available from: https://link.springer.com/article/10.1007/BF00388976. 


Marshall NB. Aspects of deep sea biology. London: Hutchinson; 1954. p. 380. 


Suga H, Ruiz-Trillo |. Development of ichthyosporeans sheds light on the origin of metazoan multicellu- 
larity. Dev Biol. 2013. https://doi.org/10.1016/j.ydbio.2013.01.009 PMID: 23333946 


Marshall WL, Celio G, McLaughlin DJ, Berbee ML. Multiple Isolations of a Culturable, Motile Ichthyos- 
porean (Mesomycetozoa, Opisthokonta), Creolimax fragrantissima n. gen., n. sp., from Marine Inverte- 
brate Digestive Tracts. Protist. 2008. 


Allen WE. Problems of Flotation and Deposition of Marine Plankton Diatoms. Trans Am Microsc Soc. 
1932 Jan; 51(1):1. 


Waite A, Fisher A, Thompson P, Harrison P. Sinking rate versus cell volume relationships illuminate 
sinking rate control mechanisms in marine diatoms. Mar Ecol Prog Ser [Internet]. 1997 Oct 16 [cited 
2021 Apr 13]; 157:97—108. Available from: http://www. int-res.com/abstracts/meps/v157/p97-108/. 


Hassett BT, Lopez JA, Gradinger R. Two New Species of Marine Saprotrophic Sphaeroformids in the 
Mesomycetozoea Isolated From the Sub-Arctic Bering Sea. Protist. 2015. https://doi.org/10.1016/j. 
protis.2015.04.004 PMID: 26046621 


Jestensen JP, Sperstad S, Johansen S, Landfald B. Molecular-phylogenetic, structural and biochemical 
features of a cold-adapted, marine ichthyosporean near the animal-fungal divergence, described from 
in vitro cultures. Eur J Protistol. 2002. 


Marshall WL, Berbee ML. Comparative morphology and genealogical delimitation of cryptic species of 
sympatric isolates of sphaeroforma (Ichthyosporea, Opisthokonta). Protist. 2013. https://doi.org/10. 
1016/j.protis.2012.12.002 PMID: 23352078 


Ducluzeau A-L, Tyson JR, Collins RE, Snutch TP, Hassett BT. Genome Sequencing of Sub-Arctic 
Mesomycetozoean Sphaeroforma sirkka Strain B5, Performed with the Oxford Nanopore minlON and 
Illumina HiSeq Systems. Microbiol Resour Announc. 2018. https://doi.org/10.1128/MRA.00848-18 
PMID: 30533723 





PLOS Biology | https://doi.org/10.1371/journal.pbio.3001551 March 29, 2022 28/31 


PLOS BIOLOGY 


Regulation of sedimentation rate shapes the evolution of multicellularity in a close relative of animals 





53. 


54. 


55. 


56. 


57. 


58. 


59. 


60. 


61. 


62. 


63. 


64. 


65. 


66. 


67. 


68. 


69. 


70. 


71. 


72. 


73. 


Grau-Bové X, Torruella G, Donachie S, Suga H, Leonard G, Richards TA, et al. Dynamics of genomic 
innovation in the unicellular ancestry of animals. elife. 2017. https://doi.org/10.7554/eLife.26036 PMID: 
28726632 


Frankel EB, Shankar R, Moresco JJ, Yates JR, Volkmann N, Audhya A. Ist1 regulates ESCRT-III 
assembly and function during multivesicular endosome biogenesis in Caenorhabditis elegans embryos. 
Nat Commun. 2017. https://doi.org/10.1038/s41467-017-01636-8 PMID: 29129923 


Dimaano C, Jones CB, Hanono A, Curtiss M, Babst M. Ist1 regulates Vps4 localization and assembly. 
Mol Biol Cell. 2008. htips://doi.org/10.1091/mbc.e07-08-0747 PMID: 18032582 


Agromayor M, Carlton JG, Phelan JP, Matthews DR, Carlin LM, Ameer-Beg S, et al. Essential role of 
hISTI in cytokinesis. Mol Biol Cell. 2009. 


Xiao J, Chen XW, Davies BA, Saltiel AR, Katzmann DJ, Xu Z. Structural basis of Ist1 function and Ist1- 
Did2 interaction in the multivesicular body pathway and cytokinesis. Mol Biol Cell. 2009. https://doi.org/ 
10.1091/mbc.e09-05-0403 PMID: 19477918 


Dudin O, Merlini L, Bendezd FO, Groux R, Vincenzetti V, Martin SG. A systematic screen for morpho- 
logical abnormalities during fission yeast sexual reproduction identifies a mechanism of actin aster for- 
mation for cell fusion. PLoS Genet. 2017; 13 (4). https://doi.org/10.1371/journal.pgen.1006721 PMID: 
28410370 


Hubner S, Jans DA, Xiao CY, John AP, Drenckhahn D. Signal- and importin-dependent nuclear target- 
ing of the kidney anion exchanger 1-binding protein kanadaptin. Biochem J. 2002. https://doi.org/10. 
1042/0264-6021:3610287 PMID: 11772400 


Hubner S, Bahr C, G6Bmann H, Efthymiadis A, Drenckhahn D. Mitochondrial and nuclear localization 
of kanadaptin. Eur J Cell Biol. 2003. https://doi.org/10.1078/017 1-9335-00308 PMID: 12800979 


Fitzky BU, Witsch-Baumgartner M, Erdel M, Lee JN, Paik YK, Glossmann H, et al. Mutations in the A7- 
sterol reductase gene in patients with the Smith-Lemli-Opitz syndrome. Proc Natl Acad Sci US A. 1998. 
https://doi.org/10.1073/pnas.95.14.8181 PMID: 9653161 


Prabhu AV, Luu W, Sharpe LJ, Brown AJ. Cholesterol-mediated degradation of 7-dehydrocholesterol 
reductase switches the balance from cholesterol to Vitamin D synthesis. J Biol Chem. 2016. https://doi. 
org/10.1074/jbc.M115.699546 PMID: 26887953 


Al-Zahrani KN, Baron KD, Sabourin LA. Ste20-like kinase SLK, at the crossroads: A matter of life and 
death. Cell Adh Migr. 2013. https://doi.org/10.4161/cam.22495 PMID: 23154402 


Rohlfs M, Arasada R, Batsios P, Janzen J, Schleicher M. The Ste20-like kinase SvkA of Dictyostelium 
discoideum is essential for late stages of cytokinesis. J Cell Sci. 2007. 


Cvrékova F, De Virgilio C, Manser E, Pringle JR, Nasmyth K. Ste20-like protein kinases are required for 
normal localization of cell growth and for cytokinesis in budding yeast. Genes Dev. 1995. hitps://doi.org/ 
10.1101/gad.9.15.1817 PMID: 7649470 


Wang Y, Wang R, Tang DD. Ste20-like kinase-mediated control of actin polymerization is a new mecha- 
nism for thin filament-associated regulation of airway smooth muscle contraction. Am J Respir Cell Mol 
Biol. 2020. https://doi.org/10.1165/rcmb.2019-03100C PMID: 31913659 


Fedotova AA, Bonchuk AN, Mogila VA, Georgiev PG. C2H2 zinc finger proteins: The largest but poorly 
explored family of higher eukaryotic transcription factors. Acta Nat. 2017. PMID: 28740726 


Qiao L, Zheng J, Tian Y, Zhang Q, Wang X, Chen JJ, et al. Regulator of chromatin condensation 1 abro- 
gates the G1 cell cycle checkpoint via Cdk1 in human papillomavirus E7-expressing epithelium and cer- 
vical cancer cells article. Cell Death Dis. 2018. 


Dasso M. RCC1 in the cell cycle: the regulator of chromosome condensation takes on new roles. 
Trends Biochem Sci. 1993 Mar 1; 18(3):96-101. https://doi.org/10.1016/0968-0004(93)90161-f PMID: 
8480369 


Hadjebi O, Casas-Terradellas E, Garcia-Gonzalo FR, Rosa JL. The RCC1 superfamily: From genes, to 
function, to disease. Biochimica et Biophysica Acta—Molecular. Cell Res. 2008. 


MO, HO, TN. The RCC1 protein, a regulator for the onset of chromosome condensation locates in the 
nucleus and binds to DNA. J Cell Biol [Internet]. 1989 [cited 2021 Sep 2]; 109(4 Pt 1):1389-97. Available 
from: https://pubmed.ncbi.nim.nih.gov/2677018/. 


MA, MWC, UV, JA. A yeast mutant, PRP20, altered in mRNA metabolism and maintenance of the 
nuclear structure, is defective in a gene homologous to the human gene RCC1 which is involved in the 
control of chromosome condensation. Mol Gen Genet [Internet]. 1990 Oct [cited 2021 Sep 2]; 224 
(1):72-80. Available from: https://pubmed.ncbi.nim.nih.gov/2277633/. https://doi.org/10.1007/ 
BF00259453 PMID: 2277633 


Forrester W, Stutz F, Rosbash M, Wickens M. Defects in mRNA 3’-end formation, transcription initia- 
tion, and mRNA transport associated with the yeast mutation prp20: possible coupling of mRNA pro- 
cessing and chromatin structure. Genes Dev [Internet]. 1992 Oct 1 [cited 2021 Sep 2]; 6(10):1914—26. 





PLOS Biology | https://doi.org/10.1371/journal.pbio.3001551 March 29, 2022 29/31 


PLOS BIOLOGY 


Regulation of sedimentation rate shapes the evolution of multicellularity in a close relative of animals 





74. 


75. 


76. 


77. 


78. 


79. 


80. 


81. 


82. 


83. 


84. 


85. 


86. 


87. 


88. 


89. 


90. 


91. 


92. 


93. 


94. 


Available from: http://genesdev.cshlp.org/content/6/10/1914 https://doi.org/10.1101/gad.6.10.1914 
PMID: 1398069 


Kadowaki T, Zhao Y, Tartakoff AM. A conditional yeast mutant deficient in mRNA transport from 
nucleus to cytoplasm. Proc Natl Acad Sci [Internet]. 1992 Mar 15 [cited 2021 Sep 2]; 89(6):2312-6. 
Available from: https://www.pnas.org/content/89/6/231 2. https://doi.org/10.1073/pnas.89.6.2312 PMID: 
1372441 


Konopka CA, Schleede JB, Skop AR, Bednarek SY. Dynamin and cytokinesis. Traffic. 2006. https://doi. 
org/10.1111/j.1600-0854.2006.00385.x PMID: 16497219 


Masud Rana AYKM, Tsujioka M, Miyagishima S, Ueda M, Yumura S. Dynamin contributes to cytokine- 
sis by stabilizing actin filaments in the contractile ring. Genes Cells. 2013. https://doi.org/10.1111/gtc. 
12060 PMID: 23679940 


Rikhy R, Mavrakis M, Lippincott-Schwartz J. Dynamin regulates metaphase furrow formation and 
plasma membrane compartmentalization in the syncytial Drosophila embryo. Biol Open. 2015. https:// 
doi.org/10.1242/bi0.20149936 PMID: 25661871 


Wang MC, Lu Y, Baldock C. Fibrillin Microfibrils: A Key Role for the Interbead Region in Elasticity. J Mol 
Biol. 2009. https://doi.org/10.1016/j.jmb.2009.02.062 PMID: 19268673 


Yin W, Kim HT, Wang SP, Gunawan F, Li R, Buettner C, et al. Fibrillin-2 is a key mediator of smooth 
muscle extracellular matrix homeostasis during mouse tracheal tubulogenesis. Eur Respir J. 2019. 
https://doi.org/10.1183/13993003.00840-2018 PMID: 30578393 


Zhang H, Apfelroth SD, Hu W, Davis EC, Sanguineti C, Bonadio J, et al. Structure and expression of 
fibrillin-2, a novel microfibrillar component preferentially located in elastic matrices. J Cell Biol. 1994. 
https://doi.org/10.1083/jcb.124.5.855 PMID: 8120105 


Balasubramanian MK, Mccollum D, Chang L, Wong KCY, Naqvi NI, He X, et al. Isolation and Character- 
ization of New Fission Yeast Cytokinesis Mutants. Genetics. 1998 Jul; 149(3):1265-75. https://doi.org/ 
10.1093/genetics/149.3.1265 PMID: 9649519 


Hirono M, Yoda A. Isolation and Phenotypic Characterization of ChlamydomonasMutants Defective in 
Cytokinesis. Cell Struct Funct. 1997; 22:1—5. https://doi.org/10.1247/csf.22.1 PMID: 9113383 


Gillmor CS, Roeder AHK, Sieber P, Somerville C, Lukowitz W. A Genetic Screen for Mutations Affecting 
Cell Division in the Arabidopsis thaliana Embryo Identifies Seven Loci Required for Cytokinesis. PLoS 
ONE [Internet]. 2016 Jan 8 [cited 2021 Sep 1]; 11(1):e0146492. Available from: https://journals.plos. 
org/plosone/article?id=1 0.1371 /journal.pone.0146492 PMID: 26745275 


Nanninga N. Cytokinesis in Prokaryotes and Eukaryotes: Common Principles and Different Solutions. 
Microbiol Mol Biol Rev [Internet]. 2001 Jun [cited 2021 Sep 1]; 65(2):319. Available from: /pmc/articles/ 
PMC99029/. https://doi.org/10.1128/MMBR.65.2.319-333.2001 PMID: 11381104 


Huang H, Fletcher L, Beeharry N, Daniel R, Kao G, Yen TJ, et al. Abnormal Cytokinesis after X-Irradia- 
tion in Tumor Cells that Override the G2 DNA Damage Checkpoint. Cancer Res [Internet]. 2008 May 15 
[cited 2021 Sep 1]; 68(10):3724—32. Available from: https://cancerres.aacrjournals.org/content/68/1 0/ 
3724. https://doi.org/10.1158/0008-5472.CAN-08-0479 PMID: 18483255 


Pfeifer F. Haloarchaea and the formation of gas vesicles. Life. 2015. https://doi.org/10.3390/ 
life5010385 PMID: 25648404 


Strand E, Jorgensen C, Huse G. Modelling buoyancy regulation in fishes with swimbladders: Bioener- 
getics and behaviour. Ecol Model. 2005. 


Villareal TA, Carpenter EJ. Buoyancy regulation and the potential for vertical migration in the oceanic 
cyanobacterium Trichodesmium. Microb Ecol. 2003. https://doi.org/10.1007/s00248-002-1012-5 PMID: 
12481233 


Sundby S, Kristiansen T. The principles of buoyancy in Marine Fish Eggs and their vertical distributions 
across the World Oceans. PLoS ONE. 2015. https://doi.org/10.1371/journal.pone.0138821 PMID: 
26465149 


Chen M, Tian LL, Ren CY, Xu CY, Wang YY, Li L. Extracellular polysaccharide synthesis in a bloom- 
forming strain of Microcystis aeruginosa: implications for colonization and buoyancy. Sci Rep. 2019. 


Millero FJ, Huang F. The density of seawater as a function of salinity (5 to 70 g kg -1) and temperature 
(273.15 to 363.15 K). Ocean Sci. 2009; 5 (2):91—100. 


Fournier ro. Observations of particulate organic carbon in the mediterranean sea and their relevance to 
the deep-living coccolithophorid cyclococcolithus fragilis [Internet]. Limnol Oceanogr. John Wiley & 
Sons, Ltd; 1968 [cited 2021 May 6]; 13:693-7. Available from: https://aslopubs. onlinelibrary.wiley.com/ 
doi/full/10.431 9/lo.1968.13.4.0693. 


Schneider C, Rasband W, Eliceiri K. NIH Image to ImageJ: 25 years of image analysis. Nat Methods. 
2012; 9:671-675. https://doi.org/10.1038/nmeth.2089 PMID: 22930834 


Wickham H. ggplot2: Elegant Graphics for Data Analysis. J R Stat Soc A Stat Soc. 2016. 





PLOS Biology | https://doi.org/10.1371/journal.pbio.3001551 March 29, 2022 30/31 


PLOS BIOLOGY 


Regulation of sedimentation rate shapes the evolution of multicellularity in a close relative of animals 





95. 


96. 


97. 
98. 


99. 


Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinfor- 
matics. 2014 Aug 1; 30(15):2114—2120. https://doi.org/10.1093/bioinformatics/btu170 PMID: 24695404 


Andrews S. FastQC—A quality control tool for high throughput sequence data. Available from: http:// 
www.bioinformatics.babraham.ac.uk/projects/fastqc/. Babraham Bioinforma; 2010. 


Smit A, Hubley R, Grenn P. RepeatMasker Open-4.0. RepeatMasker Open. 2015. 


Bao W, Kojima KK, Kohany O. Repbase Update, a database of repetitive elements in eukaryotic 
genomes. Mob DNA. 2015. https://doi.org/10.1186/s13100-015-0041-9 PMID: 26045719 


Deatherage DE, Barrick JE. Identification of mutations in laboratory-evolved microbes from next-gener- 
ation sequencing data using breseq. Methods Mol Biol. 2014. https://doi.org/10.1007/978-1-4939-0554- 
6_12 PMID: 24838886 





PLOS Biology | https://doi.org/10.1371/journal.pbio.3001551 March 29, 2022 31/31 


